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ABSTRACT 



We give a formulation of the radiative transfer equation for Lyman a photons which allows us to include the two-photon corrections 
for the 3s- Is and 3d- Is decay channels during cosmological hydrogen recombination. We use this equation to compute the corrections 
to the Sobolev escape probability for Lyman a photons during hydrogen recombination, which then allow us to calculate the changes 
in the free electron fraction and CMB temperature and polarization power spectra. We show that the effective escape probability 
changes by AP/P ~ +11% at z ~ 1400 in comparison with the one obtained using the Sobolev approximation. This speeds up of 

hydrogen recombination by AN C /N C 1.6% at z ~ 1190, implying |AC//C/| ~ 1% - 3% at / > 1500 with shifts in the positions of 

the maxima and minima in the CMB power spectra. These corrections will be important for the analysis of future CMB data. 
The total correction is the result of the superposition of three independent processes, related to (i) time-dependent aspects of the 
problem, (ii) corrections due to quantum mechanical deviations in the shape of the emission and absorption profiles in the vicinity 
of the Lyman a line from the normal Lorentzian, and (iii) a thermodynamic correction factor, which occurs to be very important. All 
these corrections are neglected in the Sobolev-approximation, but they are important in the context of future CMB observations. All 
three can be naturally obtained in the two-photon formulation of the Lyman a absorption process. However, the corrections (i) and 
(iii) can also be deduced in the normal '1 + 1' photon language, without necessarily going to the two-photon picture. Therefore only 
(ii) is really related to the quantum mechanical aspects of the two-photon process. We show here that (i) and (iii) lead to the largest 
separate contributions to the result, however they partially cancel each other close to z ~ 1100. At z ~ 1 100 the modification due to 
the shape of the line profile contributes about AN C /N C 0.4%, while the sum of the other two contributions gives AN S /N C 0.9%. 
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1. Introduction 



cosmic microwave background - early Universe — cosmology: theory - atomic processes - cosmo- 



After the seminal works of IZeldovich et alJ (Il968l) and iPeeblesl 
dl968l) on cosmological recombination, and the improvements 
in the theoretical modeling of this epoc h introduced later (e.g. 
iJones & Wyselll985t ISeager et al.ll2000h . le ading to the widely 
used standard recombination code Recfast dSeager et al.ll 1 999b . 
over the past few years the detailed physics of cosmological re- 
combination has again been reconsidered by several independent 
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groups (e.g. Dubrovich & Grachev 2005; Chluba & Sunyaev 



2006b; Khol upenko & Ivanchikl 120061; iRubino-Martm et al 



Switzer & H irata 2008; Wong & Sc ott 2007). It is clear 



that understanding the c osmological ionization history at the 
level of ~ 0.1% (e.g. see ISunvaev & Chlubal 120081; iFendt etaT] 
120081 for a more detailed overview of the different previously 
neglected physical processes that are important at this level of 
accuracy) will be very important for accurate theoretical predic- 
tions of the Cosmic Microwave Background (C MB) temperature 
and polarization angular fluctuations (e.g. see iHu etaljfl995t 
ISeliak et alj2003l) in the context of the Planck Survey otQ, which 
will be launched later this year. 

Also for a precise calibration of the acoustic horizon at re- 
combination and the possibility to constrain dark energy using 
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baryonic acoustic oscillation (e.g. [Eis£nstejnj|2005|), it is crucial 
to understand the physics of cosmological recombination at a 
high level of accuracy. Ignoring percent-level corrections to the 
ionization history at last scattering (z ~ 1100) may therefore 
also result in significant biases to the cosmolo gical parameters 
deduced using larg e catalogs of galaxies (e.g. lEisenstein et all 
120051: iHutsil 120061) as for example recently demonstrated for 
more spe culative additions to the cosmological recombination 
scenario dde Bernardis et al]|2009l) related to the possibility of 
delayed recombination dPeebles et al.ll2000l) . 

Among all the additional physical mechanisms during cos- 
mological recombination that have been addressed so far, the 
problems connected with the radiative transfer of H i Lyman a 
photons, including partial frequency redistribution and atomic 
recoil due to multiple resonance scattering, electron scattering, 
and corrections due to two-photon processes during H i recombi- 
nation (z ~ 800 - 1600), have still not been solved at full depth. 
Here we will focus on the inclusion of two-photon corrections to 
the 3s- Is and 3d- Is emission and absorption process. 

The potential importance of two-photon transitions from 
highly excited levels in hydrogen a nd helium was first pointed 
out bv lDubrovich & Grachev! d2005l) . They predicted a ~ 5% de- 
crease in the free electron fraction at z ~ 1200. However, in their 
computations of the effective two-photon decay rates for the ns 
and nd-levels they only included the first non-resonant term (i.e. 
due to the dipole matrix element connecting ns/nd — > np) into 
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the infinite sum over intermediate states. Also in their approach 
they neglected any possible transfer or reabsorption of photons 
in the vicinity of the Lyman a resonance, but simply assumed 
that all the photons accounted for by the inferred effective two- 
photon decay rate can directly escape. 

Using rate coefficients for the vacuum two-photon de- 
cay s_of_the_3s and 3 d-levels in hydrogen , as computed 



ny arogen , 

bv ICresser et al.1 d 19861), I Wong & Scottl d2007l) concluded that 
iDubrovich & Grachev! ( 2005b overestimated the impact of two- 



photon transitions on the ionization hist ory by about one ord er of 
magnitude. However, the calculation of lCresser et alJd 19861) was 
incomplete, since in their attempt to separate the '1 + 1' photon 
contributions to the two-photon formula^ from the 'pure' two- 
photon decay terms, with out clear justification they n eglected the 
first non-resonant term dChluba & Sunyaev! l2008cl) . Physically 
it seems very difficult to separate the 'pure' two-photon decay 
rate from the '1 + 1 ' resonant contribu t ions (see discussions in 
Chluba & Sunvae\fl2008ctlHiratal l2008; Ka rshenboim & Ivanovl 



20081) . e.g. because of non-classical interference effects. In a 
complete analysis these contributions s hould be taken into ac- 
count. In addition JWong & Scottl (120071) also neglected radiative 
transfer aspects of the problem. 

Slightly later this probl em was reinvestigated in more detail 
dChluba & Sunvaevl2008cl) . showing that due to two-photon de- 
cays during hydrogen recombination a decrease of more than 
~ 0.3 - 0.5% in the free electron fraction at z ~ 1150 can still 
be expected. This estimate was obtained by taking into account 
departures of the full ns-ls and nd-ls two-photon line profiles 
from the Lorentzian shape in the very distant, optically thin part 
of the red wing of the Lyman a line. In these regions it can be 
assumed that all released photons can directly escape, and hence 
lead to a successful settling of the electron in the ground state. 
No radiative transfer formulation is needed to estimate this frac- 
tion of transitions, however as mentioned in their work the cor- 
rections coming from regions with significant radiative transfer 
can still be important. According to their computations, the two- 
photon decays from s-states seem to slow hydrogen recombina- 
tion down, while those from d-states speed it up. In addition it 
was shown that the slight net acceleration of hydrogen recom- 
bination seems to be domi nated by the 3s and 3d contribution 
dChluba & Sunvaevll2008cl ). 

Another investigation of the two-photon asp ects of the re - 
combination problem was recently performed by lHiratal (120081) . 
He gave a formulation of the photon transfer problem simultane- 
ously including all two-photon corrections during hydrogen re- 
combination related to hs-1s, «d-ls, c-ls transitions and Raman 
scattering processes, also taking into account stimulated pro- 
cesses in the ambient CMB blackbody radiation field. In order to 
solve this complicated problem two approaches were used. In the 
first the two-photon continuum was discretized and turned into 
an effective multilevel-atom with virtual states related to the en- 
ergy of the photons. In the second approach the corrections were 
analytically modeled as effective modifications of the Lyman 
a and Lyman /3 decay rates. In addition, in both approaches a 
distinction between regions with '1 + 1' photon contributions 
and those with pure two-photon contributions was intr oduced 
to avo id the double-counting problem (see Sec. III.C of iHiratal 
d2008l) ) for the decay rates. As pointed out this distinction is not 
unique, but the results were s hown to be independent of the cho- 
sen parameters dHiratal2008l) . in total yielding AN e /N e ~ +13% 
at z ~ 900 and ANJN e ~ -1.3% at z ~ 1300. 



Given the delicate complexity of the two-photon transfer 
problem it is very important to independently cross-validate the 
results obtained by different groups. In this paper we offer an- 
other approach to this problem in which we take into account 
the two-photon nature of the 3s- Is and 3d- Is decay channels, 
without introducing any criterion distinguishing between 'pure' 
two-photon decays and '1 + 1' resonant contributions. We give 
a formulation of modified rate equations for the different hydro- 
gen levels and the photon transfer equation, which we then use 
to compute the effective H i Lyman a photon escape probability 
including these corrections. 

Although it is clear that in particular the atomic recoil 
effect speeds up hydrogen recombination up at the percent- 
level dGrachev & Dubrovichl 120081: IChluba & Sunvaevfl2008bl) 
and also partial frequency redistribu tion will lead to some ad- 
ditional modification^ here like in IHiratal d2008l) we will ne- 
glect the frequency redistribution of photons due to resonance 
scattering and work in the no line sc attering approximation. As 
explained in several previous works ( Chluba & Sunyaev 2008b; 
ISwitzer & Hiratall2008l:lRubino-Martm et al. 2008; Hirata 2008b 
for conditions in our Universe (practically no collisions) this is a 
much better description than the assumption of complete redis- 
tribution, which is used in the derivation of the Sobolev escape 
probability. We also take into account stimulated 3s- Is and 3d- 
ls two-photon emission, finding this process to be sub-dominant. 
However, until now we do not include the effect connected with 
Raman scattering in this paper. 

Instead of solving the obtained coupled system of equations 
simultaneously, we assume that the corrections will be small, so 
that each of them can be considered as a perturbation of the nor- 
mal '1 + 1' photon result. Therefore we can use pre-computed 
solution^ for the populations of the different hydrogen levels 
as a function of time to obtain the time-dependent photon emis- 
sion rate for the different decay channels. This approach allows 
us to solve the Hi Lyman a radiative transfer equation semi- 
analytically also including the 3s- Is and 3d- Is two-photon cor- 
rections. Using the obtained solution for the spectral distortion 
at different redshifts one can then compute the effective Lyman 
a escape probability as a function of time. This value can be di- 
rectly compared to the normal Sobolev escape probability which 
then also allows to deduce the expected modification in the cos- 
mological ionization history and CMB temperature and polar- 
ization power spectra. 

Here we show that the effective escape probability changes 
by AP/P — 1-11% at z ~ 1400 in comparison with the value de- 
rived in the Sobolev approximation (see Fig. |T7t. As we explain 
in detail, this total correction is the result of the superposition 
of three independent processes, related to (i) time-dependent as- 
pects of the problem, (ii) corrections due to deviations in the 
shape of the emission and absorption profiles in the vicinity of 
the Lyman a line from the normal Lorentzian, and (iii) a thermo- 
dynamic correction factor. All these corrections are neglected in 
the cosmological recombination problem, but for the analysis of 
future CMB data they should be taken into account. 



This expression was first derived bv Goppert-M aveJ j 193 lh 



3 As already mentioned in IChluba & Sunvaevl d2008bl) . our current 
version of the code already contains the corrections due to line diffu- 
sion on therma l atoms, atomic recoil and ele ctron scattering. In good 
agreement with Grachev & Dubrovichl d2008h we found that atomic re- 
coil is most important, but partial freq uency redistribution only l eads 
to an additional small modification. In IChluba & Sunvae\J d2009h we 
will present the results of these computations, also taking the 3d- Is and 
3s-ld two-photon corrections into account. 

4 We used the output of the latest version of our multilevel code 
jRubino-Martfn et alj|2006l : IChluba et all2007l) . 
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In the '1 + 1' photon picture the pur ely time-dependent 
correct ion was already discussed earlier (IChluba & Sunvaevl 
I2008bl) . showing that changes in the state of the medium (e.g. 
number densities and Hubble expansion rate) cannot be ne- 
glected in the computation of the Lyman a escape probability. 
This is because only a very small fraction (~ 10~ 4 - 10~ 3 ) of 
all interactions with the Lyman a resonance lead to a complete 
redistribution of photons over the whole line profile. As a conse- 
quence only the region inside the Doppler core reaches full equi- 
librium with the photon occupation number at the line center and 
can be considered using quasi-stationary conditions. However, 
outside the Doppler core time -dependent aspects of the p roblem 
have to be taken into account (IChluba & Sunvaevll2008bl) . 

The second correction is related to quantum mechanical 
modifications in the shape of the line profiles describing the ns- 
ls and «d-ls two-photon decay channels. As we explain here, 
this is the only correction that can only be obtained when us- 
ing the two-photon pictur e. As already discussed earlier (e.g. 
IChluba & S unyaev 2 008d) . this leads to deviations of the corre- 
sponding profiles from the normal Lorentzian. One consequence 
of this is that, depending on the considered process, more (for 
«d-ls transitions) or fewer (for ns-ls transitions) photons will 
directly reach the very distant red wing (xq < -1000 Doppler 
width), where they can i mmediately escape. This cor rection was 
already estimated earlier dChluba & Sun yaev 2008c), but here it 
will now be possible to refine these computations, also extend- 
ing to regions closer to the line center, where radiative transfer 
effects are important. Similarly, modifications in the blue wing 
emission can be taken into account using the approach presented 
here. Most importantly, because of the correct inclusion of en- 
ergy conservation, the two-photon profiles will not extend to ar- 
bitrarily high frequencies. This will avoid the low redshift self- 
feedback that was recently seen in a time-dependent formul ation 
of the Lyman a escape problem dChluba & Sunvaevl2008bl) . and 
which here can be modeled more consistently. 

The last and also most important correction discussed in this 
paper is related to a frequency-dependent asymmetry between 
the line emission and absorption process, that is normally ne- 
glected in the deriv ation of the Sobolev escap e probability. As 
pointed out earlier dChluba & Sun yaev 2008b), within the nor- 
mal ' 1 + 1 ' photon formulation for the line emission and absorp- 
tion process especially in the damping wings of the Lyman a line 
a blackbody spectrum is not exactly conserved in full thermody- 
namic equilibrium. This leads to the requirement of an additional 
factor, f v , inside the absorption coefficient, which in the '1 + 1' 
photon picture can be deduced using the detailed balance prin- 
ciple (see AppendixlBli. However, within the two-photon formu- 
lation this correction naturally appears in connection with the 
two-photon absorption process, where one photon is taken from 
close to the Lyman a resonance and the other is drawn from 
the ambient CMB blackbody photon field at frequency follow- 
ing from energy conservation^ (see Sect. 12.1. H and in particular 
Sect. l33^1 l. 

We will henceforth refer to f v as the thermodynamic cor- 
rection factor. It results in a suppression of the line absorption 
probability in the red, and an enhancement in the blue wing of 
the Lyman a resonance. This asymmetry becomes exponentially 
strong at large distances from the resonance. In most astrophys- 
ical applications one is not interested in the photon distribution 



5 For the ls-3d two-photon absorption this will be V = V31 - v, where 
V3i is the corresponding ls-3d transition frequency and v denotes the 
frequency of the photon that is absorbed in the vicinity of the Lyman a 
resonance. 



very far away from the Lyman a line center, so that this correc- 
tion usually can be neglected. However, for the cosmological re- 
combination problem even details at distances of ~ 1 3 - 10 4 
Doppler width do matter dChluba & Sunvaevl l2008bl) . so that 
such an inconsistency in the formulation of the transfer problem 
has to be resolved. As we will show here the associated correc- 
tion is very important, leading to a significant speed-up of H 1 
recombination. 

We also demonstrate that including all three modifications 
to the escape probability, the number density of free electrons 
is expected to change by AN e /N e ~ -1.3% (see Fig. [18). 
close to the maximum of the Thomson visibility function 
dSunvaev & Zeldovich|[T970h at z ~ 1100, which matters most 
in connection with the CMB power spectra. The 3s- Is and 3d- 
ls two-photon corrections (related to the shape of the profiles 

and the thermodynamic factor alone) yield AN e /N e 2.4% at 

Z ~ 1110. A large part (~ 1.1% at z - 1100) of this correction 
is canceled by the contributions from the time-dependent aspect 
of the problem (see Fig. [T8l for detail s). Our results seem to be 
rather similar to those of lHiratal (2008) for the contributions from 
high level two-photon decays aloneQ 

We also compute the final changes in the CMB temperature 
and polarization power spectra when simultaneously including 
all processes under discussion here (see Fig. |T9l >. The correc- 
tions in the E-mode power spectrum are particularly impressive, 
reaching a peak to peak amplitude of ~ 2% - 3% at / > 1500, 
and significant shifts in the positions of the maxima in the CMB 
power spectra. Taking these corrections into account will be im- 
portant for the future analysis of CMB data. 

The paper is structured as follows: in Sect. [2] we give the 
equation for the modified Lyman a transfer problem. There we 
infer the equations by generalizing the normal '1 + 1' photon 
transfer equation in order to account for the mentioned pro- 
cesses. In the Appendix lAl we give a more rigorous derivation 
using the two-photon formulae, also generalizing the rate equa- 
tions for the different hydrogen levels. We then give the solution 
of the transfer equation in Sect. 12.21 and show how to use it to 
compute the effective Lyman a escape probability (Sect. 12.31 ). 
We explain the main physical differences and expectations for 
the corrections in comparison with the '1 + 1' photon formu- 
lation in Sect. [3] Then we include 'step by step' the different 
correction terms and explain the changes in the results for the 
spectral distortion around the Lyman a line (Sect. [4) and the ef- 
fective escape probability (Sect. [5]). In Sect. [6] we then give the 
results for the ionization history and the CMB temperature and 
polarization power spectra. We conclude in Sect. [7] 



2. Two-Photon corrections to the Lyman a emission 
and absorption process 

The aim of this Section is to write down the line-emission and 
absorption terms describing the evolution of the photon field 
in the vicinity of the Lyman a resonance including the 3s- Is 
and 3d- Is two-photon corrections. Here we will try to motivate 
the form of this equation in terms of the additional physical as- 
pects of the problem that should be incorporated. We refer the 
interested reader to Appendix [A] in which we provide the ac- 
tual derivation of this equation using a two-photon formulation. 
There the central ingredient is that the photon distribution around 
the Balmer a line is given by the CMB blackbody. This fact 



6 Note that this is o nly part of the total correction which was pre- 
sented in lHiratj j2008h . 
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makes it possible to rewrite the two-photon transfer equation as 
effective equation for one photon, as presented here. 

In this Section we also give the solution of the modified 
transfer equation (Sect. 12.21 ) and explain how one can use it to 
compute the effective escape probability for the Lyman a pho- 
tons (Sect. |23). 



the photon distribution in the vicinity of the second photon that 
is involved in the Lyman a absorption process (Sect. [33^21 . This 
is due to the fact that the photon which enables the 2p-electron 
to reach the 3s, 3d, or continuum is drawn from the ambient ra- 
diation field, which in the cosmological recombination problem 
is given by the CMB blackbody. 



2.1. Modified equation describing the emission and death of 
Lyman a photons 

Within the semi-classical formulation of the Lyman a transfer 
equation every relevant physical process is envisioned as a sin- 
gle step process involving one photon of the photon field. This 
leads to the introduction of photon death an d scattering proba- 
bilities that only depend on redshift (e.g. see lChluba & Sunvaevl 
l2008bl) . Also in the single photon picture the line profiles for 
the different Lyman a emission and absorption channels under 
the assumption of complete redistribution are all identical. For 
example, it will not make difference if the electron reaches the 
2p-state and then goes to the 3s, 3d or continuum. In all three 
cases the absorption profile will be given by the usua l Voigt- 
profile. As explained earlier dChluba & Sunvaevll2008bl) . in the 
normal '1 + 1' photon language the Lyman a line-emission and 
absorption terms can be cast into the form 



I ^1 
c dt 



Ly-tr 



An Av d 



p/ kV21 Bn 



N ls N v ]. (1) 



Here (py(v) is the usual Voigt-profile (see Appendix iDl for def- 
inition), with normalization JT°° dvdO = 1, where Av D is 



the Doppler- width of the Lyman a line. Furthermore, pl^ is the 
Lyman a emission probability in the ' 1 + 1 ' photon pictur^H as 
given by Eq. (|A.7t , and p d y — 1 - pl^ the corresponding death 
probability. R^ p describes the rate at which fresh electrons are 
added to the 2p-state, and is defined by Eq. ( |A.9bt . 



2.1 .1 . Introducing the thermodynamic correction factor 

As mentioned in the introduction, in this form Eq. fl} does not 
exactly conserve a blackbody spectrum in the case of full ther- 
modynamic equilibrium. Knowing the ' 1 + 1 ' photon line emis- 
sion term and using the detailed balance principle one can obtain 
the thermodynamic correction facto^ 



Mz) = 



h[v-v 21 ]/kTJ;) 



(2) 



which is n ecessary to avoid this prob lem. This factor was in- 
troduced in lChluba & Sunvaevl d2008bl) already. Inserting it into 
Eq. ([1]) we then have 



c dt 



Ly-a 



An Av D 



[Pern R 2p 



p l d 7 hV2 



BnNufyNy]. (3) 



In the standard ' 1 + 1 ' photon formulation f v has no direct phys- 
ical interpretation. It is simply a consequence of thermodynamic 
requirements on the form of the equations. However, as shown in 
Appendix|A]the same factor f v naturally appears in a two-photon 
formulation of the problem. It is actually related to the shape of 



7 For formal consistency we included the factors 1 + n p \(vii) due to 
stimulated emission in the ambient CMB blackbody field in the defini- 
tion of p e £, although during recombination /t p i(v2i) * 1- 

8 We added a short derivation for f v in Appendix iBl 



2.1 .2. Including the corrections due to the profiles of the 
different decay channels 

As a next step we want to take the differences between the line 
profiles of the different absorption and emission channels into 
account. One can see that in Eq. ([3]) there is no distinction made 
between the different routes the electron took before or after en- 
tering the Is <-> 2p transition. However, as mentioned in the 
introduction, the line-emission profiles depend on how the fresh 
electron reached the 2p-state via channels other than the Lyman 
a transition. 

In order to distinguish between the different possibilities 
(e.g. Is «-» 2p <-> 3s/3d/c), one should allow for profiles, 0,(v), 
that depend on the channel ;. Also the partial rate at which elec- 
trons enter the 2p-state will depend on i, leading to the replace- 
ment Ri — > R': + with Rt = Y: R': + , where the sum runs over all 

zp -p zp 

possible ' 1 + 1' photon channels via which the number of Lyman 
a photons can be affected. Furthermore, the probability with 
which electrons are absorbed will become channel-dependent, 
so that p 1 / -> p' d with pi 7 = £ ; p\. 

Here it is important that F£j* and p' d both will only depend on 
time but not on frequency. This is because microscopically it is 
assumed that the absorption process leads to a complete redistri- 
bution over the profile 0,(v). With this comment it is also clear 
that the factor f v should be independent of the channel, since 
otherwise detailed balance for each process cannot be achieved. 

With this in mind it is clear that the more general form of 
Eq. (O should read 



1 



cFV, 
dt 



Tj rr t& 4 + - /1V21 B 12 /» N A ■ (4) 



Ly 



An Av D 



In AppendixlAlwe argue that both R'£ and p' d can be given using 
the normal '1 + 1' photon values for the different rates. We also 
specify how to compute the profiles, <pi(v), including stimulated 
two-photon emission (Sect. [Cj. However, in what follows it is 
only important that non of these depends on the solution of the 
problem for the photon field. This is because we assume that the 
readjustments in the populations of the different level or number 
density of free electrons is small and hence can be neglected to 
lowest order. Numerically one can include the correction to the 
correction iteratively, but we leave this for a future paper. 

It is important to mention that because for two-photon tran- 
sitions ns/d — > Is from « > 3 also photons connected with 
the other Lyman series are emitted, Eq. © in principle can be 
used to describe the simultaneous evolution of all Lyman se- 
ries photons. Similarly, one can account for the two-photon cor- 
rections due to transitions from the continuum c — > Is, simul- 
taneously including the Lyman continuum and all other con- 
tinuua. However, in this case one can no longer clearly dis- 
tinguish between the different Lyman series. Also the equa- 
tion will simultaneously d escribe the process of Ly-n feedback 
dChluba & Sunvaevll2007l) . in addition accounting for its exact 
time-dependence. To avoid these complications, below we will 
first only take into account the two-photon corrections for the 
3s- Is and 3d- Is channel, but leave the others unchanged. In this 
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over the absorption profile 

z) 



case it is possible to directly compare the results with the Lyman 
a problem. In Sect. [7] we briefly discuss the expected effect of 
this approximation, but leave a detailed analysis for another pa- 
per. 



2.2. Solution of the transfer equation 

For a given ionization history, the solution of Eq. <j4j in the ex- 
panding Unverse can be readi l y found , using the proceedure de- 
scribed in lChluba & Sunvaevl ( l2008bl) . If we introduce the effec- 
tive absorption optical depth as 



T a bs( 



v,z',z) = j 



Pd 



ccr r N 



H(l+Z) 



bs (x[l+z],z)dz (5a) 



s (v, z) = f v (z) V Mv, z) = f v (z) 4>l bs (v, z) (5b) 
i Pd 



with - lr 



^r 1 it 1 - and the dimensionless fre- 

4tt Avd 



P d ' = 1 - Pem, CT, 

quency x = v/(l + z), and define the effective emission profile 



<Pem(v, Z) = ^ — f </>i(V, Z) 



2p 



then Eq. © takes the simple form 

= /#V r N u abs (v, z) {N e v m - N v ] 



I ^1 

c dt 



Ly 

em_ 2v 21 gH 



X . . — *»em 



C 2 g2pR 2v Nls 0abs(v,z) 



0abs(v, z) ' 



(6) 

(7a) 
(7b) 



where N em is only redshift dependent. 

The solution of t his equation in the expandin g Universe can 
be directly given (see lChluba & Sunvaevll2008bl) 



AN? ym (z) = [N em (z) - <] x F(v,z). 



(8a) 



Here the function F v represents the frequency dependent part of 
the solution for the spectral distortion, which is defined by 



F(v,z) =j& 



(z,z')d z >e 



-T abs (v,:',z) 



dz' 



© a (z,z') 



N em (z') x 



n(vV) 



„(W) 



N em (z) - N? 21 
, NUz')x^m-N p ! 



(*ab s ( V '^') 



fy 



N em (z) - Nt 



(8b) 
(8c) 

(8d) 



where AN V = N v - N*\ V = jc[1 + z'] and at z > z s the CMB 
spectrum is assumed to be given by a pure blackbody spectrum 
N?. Furthermore, A> m (z) = N em (z)/[1 + z] 2 , x 21 = y 21 /[l + z], 
4i = v 2 i/[l +z'], and Jvf = ^je^j, with T = 2.725 K. Note 
that A/£ does not explicitly depend on redshift. Also we have 
used that fAz'Wf = N P J . 



2.3. Number of absorbed photons and the effective Lyman a 
escape probability 

With the solution (0 one can directly compute the number of 
absorbed photons. For this we define the mean of N v = I v /hv 



Kt , , f 0abs(v,z) Ar f ^ab S ( V ' Z ) , .r , , n 

A^abs(z) = -; — t N v dvdQ = — — f v N v dvdQ. 

J 47rAv D J 47rAv D 

= < ^ab s + / Z) /v(z) AiV v dv. 



(9) 



where we have set y?* bs (v, z) = < KbS v,z ^l ^ Vd an< ^ introduced the 
norm of y?* bs as ^* bs = J ip* bs (v,z)dv. If we now insert the solu- 
tion (0 into this expression we can write 



AN'2 m (z) = [N em (z) - <,][! - P] 



P = 1 



j </4s( V )/v 



(z)F v dv. 



(10a) 
(10b) 



Here P will later be interpreted as the main part of the effective 
escape probability (see Sects. I2~3~!2"l andl5"ll. 
Similar to iV a bs(z) one can also define 



J 4ttAv d J 4ttAv d 

— A^em ^em > 

so that with the transfer equation (0 it follows 



N em (z)dvd£l 
(11) 



dA y 
~~dT 



dvdQ. 



Ly 



~ c J dr 



= p^hv2iB n Nu[N, 



Ly 



em ^em 



-[N em -N? 2l ][l-P]} 
= P 1 J hv 2 iBn N u {P eff [A em - <j} , (12a) 



with 



P eff = P + AF ind 

A^em A^em 



APind = 



■< ^bs 



iV em iV V2i 



(12b) 
(12c) 



where Aip em = <p em - 1 and A^* bs = <io* bs - 1 . As we explain below, 
with these definitions the effective escape probability, P e ff, can 
now be directly compared with the value in the normal ' 1 + 1 ' 
photon formulation and the Sobolev escape probability. 

2.3.1 . Range of integration over the profiles 

In the above derivation we have not specified the range of inte- 
gration. Since the 3s and 3d two-photon profile include both the 
Balmer a and Lyman a photons, by carrying out the integrals 
over the frequency interval (0, oo), one would count 2y per tran- 
sition. In order to avoid this problem, we can simply restrict the 
range of integration to v > vn/2, but leave all the other defini- 
tions unaltered. Since V31/2 is far away from the Lyman a res- 
onance this does not lead to any significant problem regarding 
the normalization of the normal Voigt-functiorQ. In addition, for 
the quasi-stationary approximation the contribution to the value 
of the escape probability from this region are completely negli- 
gible. Therefore this restriction does not lead to any bias in the 
result, but does simplify the numerical integration significantly. 



9 V31/2 corresponds to x D « -1.7 X 10 4 [yj^j] Doppler width, so 
that the absolute error in the normalization of V is ~ 1.6 x 10~ 8 . 
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2.3.2. Relating the corrections in the spectral distortion to the 
corrections in the effective escape probability 



10-shell Hydrogen Atom 



and ^ relate 



We now want to understand how differences in F 
to corrections in the effective escape probability. For this we first 
want to emphasize that in the normal ' 1 + 1 ' photon picture, under 
the assumption of quasi-stationarity and in the no line scattering 
approximation, following the derivation of the previous Section 
one would find (IChluba & Sunvaevll2008bh 



~dT 



p'jhvnBnNuPd [N a 



(13) 



Ly 



with P d = 



and r d = p^rs 



diV d 

It is clear that At represents the effective change in the 
total number density of photons involved in the Lyman a evolu- 
tion over a short time interval At, and hence is directly related 
to the change in the total number of electrons that settle in the 
ground state via the Lyman a channel. By comparing Pd with 
P e ff, as defined by Eq. d!2bt . one can therefore deduce the re- 
quired effective correction to the Sobolev escape probability, that 
is normally used in the formula tion of the recombination p rob- 
lem. Following the arguments o flChluba & Sunvaevl ( l2008bl) this 
correction should be given by 



AP S 



P d Pett 



Ps, 



(14) 



where Ps 



is the standard Sobolev escape probability, 



with the usual Sobolev optical depth, rs. 

3. Main sources of corrections to the Lyman a 
spectral distortion 

Using the solution © one can already identify the main sources 
for the corrections to the photon distribution in comparison with 
the quasi-stationary approximation. These can be split up into 
those acting as a time and frequency dependent emissivity, which 
is characterized by a , and those just affecting the absorption 
optical depth, r a b s . Below we explain how the two-photon as- 
pect of the problem enters here, and which effects are expected. 
In Sect.[4]and[5]we discuss the corrections to the Lyman a spec- 
tral distortion and the effective escape probability in comparison 
with the standard '1 + 1' photon formulation in more detail. 



3.1. Relative importance of the different Lyman a absorption 
channels 

Before looking at the solution of the transfer equation in more 
detail it is important to understand, which channels on average 
contribute most to the absorption of Lyman a photons. In Fig.Q] 
we present the partial death probabilities for different channels, 
as defined in the Appendix lAl At all considered redshifts more 
than ~ 90% of the absorbed Lyman a photons disappear from 
the photon distribution in ls-3d two-photons transition. In con- 
trast to this, only about 2% of all transitions end in the 3s-state. 
This is because the ratio of the 2p-3s and 2p-3d transition rates is 
about g3 S A-} S 2p/g3dA-i S 2p ~ 1/50. One can also see that in general 
the ls-«d channels are more important than the ls-«s channels, 
and that the contributions of ls-3s and ls-4d two-photon chan- 
nels are comparable, where at high redshifts the ls-4d channels 



all channels 
continuum 
s-state 
d-state 

other channels 




Fig. 1. The death probabilities for different Lyman a absorp- 
tion channels. We used a 10-shell hydrogen atom. The thick 
solid line shows the total death probability, p d (for definition see 
Appendix [A). 



contributes slightly more (~ 2% versus ~ 7%). However, less 
than ~ 0.5% of photons are directly absorbed to the continuum. 

Assuming that the final modification in the ionization history 
is AN e /N e ~ 1% when only including the two-photon aspects for 
the 3d- Is channel, then the above numbers suggest that: (i) the 
additional correction is expected to be similar to AN e /N e ~ 0. 1 % 
when also taking the two-photon character of the ls-3s, ls-4d, 
and ls-5d channels into account; (ii) neglecting the two-photon 
character for the transition to the continuum should lead to an 
uncertainty of AN e /N e < 0.1%. These simple co nclusions seem 
to be in good agreement with the computations of lHiratal d2008l) . 
This also justifies the fact that here as a first step we only con- 
sider the two-photon corrections to the 3s- Is and 3d- Is channel. 
However, we plan to take the other two-photon corrections into 
account in a future paper. 

3.2. Effective Lyman a emission and absorption profile 

As we have seen in the previous section, the main channel 
for Lyman a absorption is due to the ls-3d two-photon tran- 
sition. This implies that the effective absorption profile, 0* bs , 
will be very close to the one following from the 3d- Is chan- 
nel alone. In Fig.|2]we give the spectral dependence of different 
line profiles in the vicinity of the Lyman a resonance at red- 
shift z = 1300. For comparison we also show the Voigt profile, 
(py (see Appendix iDb. One can clearly see the asymmetry of the 
two-photon profiles around the Lyman a line center and the de- 
viations from the Lorentzian shape in the distant damping wings. 

In the right panel we also show the effective emission profile, 
<p em , for the 3 shell atom, as defined by Eq ((5). In the computa- 
tions we only included the 3s and 3d two-photon profiles, but as- 
sumed that in the continuum channel ( 1 s <-> 2p <-> c) photons are 
emitted according to the normal Voigt profile. As one can see the 
effective emission profile indeed is very close to the 3d- Is two- 
photon profile, including stimulated emission. Only at v > V31 
one can see the small Lorentzian contribution from the contin- 
uum channel. Close the V31 one can also see the small admixture 
of the 3s- Is two-photon profile. As can be deduced from the left 
panel in Fig. [2] at v ~ V31 the stimulated 3s-ls two-photon pro- 



Chluba and Sunyaev: Ly a escape during cosmological hydrogen recombination 7 




Fig. 2. Different line profiles in the vicinity of the Lyman a resonance at redshift z = 1300. The left panel shows the 3s and 3d 
emission profiles in comparison with the normal Voigt profile. In the right panel we show the effective emission profile for a 3 
shell hydrogen atom, as defined by Eq. (O, in comparison with the 3d emission profile and the Voigt profile. The curves labeled 
'induced' include the effect of stimulated two-photon emission due to enhancement connected with the second photon released at 
low frequencies. This process is only important close to v ~ V3 1 and eventually leads to a sub-dominant correction of AN e /N e < 0. 1 % 
in the ionization history. On the other hand, the overall asymmetry in the e ffective emission profile (re d wing stronger that blue wing) 
has important implications for the hydrogen recombination problem (see lChluba & Sunvaevll200"8a) . 



file is about ~ 40 times larger than the 3d- Is two-photon profile. 
With appropriate renormalization one can also obtain this fac- 
tor using the approximation ( IC.3b . Although ~ R^ + /50, 

due to this factor at v ~ V31 the 3s channel adds about | 03d, or 
~ 44% to the effective emission profile. 

3.3. Time and frequency dependence of the absorption 
optical depth 

In the definition of F(v, z), Eq. d8bl . the function a accounts 
for the frequency and time dependence of the emission process. 
For a = 1 the shape of the solution for the spectral distortion 
depends only on the absorption optical depth, r a b s , as defined by 
Eq. ( l5at . In this case one can directly write 

F Q (v,z) = 1 - e - T ** ivM \ (15) 

Separating this part of the solution turns out the be very useful 
for numerical purposes. However, as we will see in Sect. 13.4.21 
Fq does not describe the main behavior of the spectral distortion 
when including the thermodynamic correction factor /,,. 

3.3.1 . Pure time-dependent correction to T a b S 

If we neglect the two-photon corrections to the 3s and 3d pro- 
files ((pi = <py) and set f v = 1 then we can look at the 
pure time-dependent corre ction to r a b s . As explained earlier 
dChluba & Sunvaevll2008bl) . the dependencies of p^, N\ s , and H 
on redshift lead to deviations of the solution for the spectral dis- 
tortion from the quasi-stationary case. Here the most important 
aspects are that, depending on the emission redshift, the total 
absorption optical depth until the time of observation (here z), is 
effectively lower (for z em > 1400), or greater (for z em S= 1400) 
than in the quasi-stationary case. In addition the deviation from 



the quasi-stationary case depends on the initial frequency of the 
considered photon, since close to the line center photons travel 
a much shorter distance before getting absorbed than in the very 
distant wings, implying that time-dependent corrections are only 
important for photo ns that are emitted outside t he Doppler core 
(for more details see lChluba & Sunvaevll2008bl) . 

In Fig. [3] we illustrate these effects on T a b s for emission red- 
shift z = 1100. We show the optical depth as a function of the 
initial frequency for different Az. In the upper panel we show the 
results for the case under discussion here (solid line). For com- 
parison we show the values of the optical depth using the normal 
quasi-stationary optical depth (dashed lines) for which one has 

r*V j ~ 

7f(v,Z,Z , )*T d (z) <Mv)-A (16) 

with Ta(z) = P^Ts, where ts is the normal Sobolev optical 
depth, V = vjftr and z' - z - Az. 

For very small Az/z one expects no significant difference be- 
tween the full numerical result for r a b s and this approximation. 
However, looking at the cases Az/z = 10~ 5 , 10~ 4 and 10~ 3 one 
can see that even then there is a small difference in the distant 
red and blue wings of the line. This is not due to time-dependent 
corrections but due to the fact that, as usual, in Eq. ( fT6] l we ne- 
glected the factor V21/V which appears in the definition of r a bs, 
leading to r a b s /Tj S ~ V21/V > 1 on the red, and Tabs/tr < 1 on 
the blue side of the resonance. 

For the cases Az/z = 0.01 and 0.1 we start to see the cor- 
rections due to the time-dependence. Here most interestingly for 
Az/z = 0.1 in both wings r a b s < r^ s , is because the photons 
were releas ed at z < 1400, so that tj(z) decreases while the pho- 
tons travel dChluba & Sunyaev 2008b). This means that the in- 
tegral over different redshifts t<j(v, z,z') ~ J yl Td(z)0v(v)^ can- 
not reach the value for t] s (z). Note that comparing with the value 
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3-shell Hydrogen Atom, : 




-5000 
X D 

3-shell Hydrogen Atom, 



10000 




-5000 
X D 

3-shell Hydrogen Atom, z 




Fig. 3. Modifications in the absorption optical depth 
T a bs(*D,Z,Z - Az) for z = 1100. Here x D = [v - V2i]/Av D , 
where Avd is the Doppler width of the Lyman a resonance. In 
each plot we show a sequence (lower to upper set of curves) of 
Az/z = 1(T 5 , 10~ 4 , 10~ 3 , 10~ 2 and 0.1. For detailed explanation 
see Sect. [3] 



Tj S (z') < Tj S (z) at the absorption redshift z' — z— Az < z one finds 
Td(v, z, z') > t^(z') following a similar argument. Usually this is 
the comparison which is made when talking about the escape 
probability at redshift z, so that the role of z and z' is simply 
interchanged. 

The difference due to the time-dependence is not yet very 
visible for Az/z = 0.01 (the changes should be \At/t\ ~ \Az/z\), 
but one can see it in the region < xu < xjj ~ 500. There 
it is clear that the emitted photons will reach the Doppler core 
over a period that is shorter than the chosen Az/z. For the case 
Az/z = 0.1 this region is < xu < ~ 4800. Depending on 
how far the photon initially was emitted from the Doppler core 
the time it will travel before reaching ~ will grow with 
increasing x^. This implies that at the redshift z c < z < 1400 of 
Doppler core crossing Td(z c ) < Td(z), leading to the slope seen 
in the regions < xd ^ x c D . 

Note that in the final result the time-dependent correction 
to Tabs is not so important, only leading to modifications in the 
escape probability by \AP/P\ ~ 1% - 3%. The time-dependence 
of 0' is much more relevant (see Sect.[5]for more details). 



3.3.2. Effect of Thermodynamic correction factor on r abs 

If we now include the thermodynamic correction factor f v , as 
given by Eq. @, in the computation of T a b s , then it is clear 
that for photons appearing at a given time on the red side of 
the Lyman a resonance, the total absorption optical depth over 
a fixed redshift interval will be lower than in the standard ap- 
proach, independent of the emission redshift. Since liV2\/kT y ~ 

40 fTTro] 1 one has " y 2i]A7V ~ p [i^] _1/2 . Due to the 
exponential dependence of on the distance to the line center 

this implies that at xd < -10 3 [yuij] ' photons even directly 
escape, without any further reabsorption. This is in stark con- 
trast to the standard approximation (f v = 1) for which even at 

distances 10 4 some small fraction of photons (comparable to 

10 -3 at z ~ 1 100) still disappears. We illustrate this behavior in 
the central panel of Fig. [3] where at large distances on the red 
side of the resonance the value of r a bs is many orders of magni- 
tude smaller than in the quasi-stationary approximation. As we 
will see below (e.g. Sect. 15. U . the thermodynamic factor leads 
to the largest correction discussed in this paper, and in fact it 
is this red wing suppression of the absorption cross section that 
contributes most. 

As mentioned in Sect. 12.1.11 physically this behavior re- 
flects the fact that the photon which enables the 2p-electron to 
reach the 3s and 3d is drawn from the ambient CMB radia- 
tion field. For photons on the red side of the Lyman a reso- 
nance (v < V21) a photon with V > V32 is necessary for a Is 
electron to reach the third shell. Since during Hi recombination 
the Balmer a line already is in the Wien tail of the CMB, this 
means the that relative to the Balmer a line center the amount 
of photons at V > V32 is exponentially smaller, depending on 
how large the detuning is. Denoting the frequency of the sec- 
ond photon (absorbed close to the Balmer a resonance) with 
V = V31 - v, by taking the ratio of the photon occupation num- 
bers ri/n pl {v3 2 ) ~ « P i(v3i - v)/n p i(y 32 ) * e hlv ' V2 ' ]/kT -> we again 
can confirm the exponential behavior of f v . Note that the same 
factor will appear even when thinking about two-photon transi- 
tions towards higher levels with « > 3 or the continuum. It is 
a result of thermodynamic requirements, which should be inde- 
pendent of the considered process, as long as the second photon 
is drawn from the CMB blackbody. 
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On the other hand, for photons that are released on the blue 
side of the Lyman a line the the total absorption optical depth 
is larger than in the standard approximation (see Fig. [3] central 
panel for illustration). Due to the exponential dependence of f v 
on frequency for 0* bs = <py this even leads to an arbitrarily large 
absorption optical depth in the very distant blue wing. Again 
this behavior can be understood when thinking about the second 
photon as drawn from the CMB blackbody. However, now there 
are exponentially more photons available than at the Balmer a 
line center. 

This very strong increase in the absorption optical 
depth implies that photons are basically reabsorbed quasi- 
instantaneously, so that the usual quasi-stationary approxima- 
tion for the computation of r a b s should be possible, like inside 
the Doppler core. In this case one therefore has 



T q Uv,Z s ,Z) * T d (z) 



, vi\ dv' 

0abs(v ) 

AVn V 



(17) 



with TdCz) = p 1 / t$, where t$ is the normal Sobolev optical 

depth, and v s = v-r^. This approximation for r a b s will also be 
very accurate close to the line center, but is very crude in the 
red wing. Note that for f v — 1 and <f>* dbs = <pv, tL recovers 
the approximation for the normal absorption optical depth in the 
quasi-stationary approximation, Eq. ( fToT i. 

For z s — > oo, y » V21 and assuming that 0* bs = <py one has 



T^ s e (v,z s ,z)*T d (z)Mz)- 



7TX D 



(18) 



With this equation it is possible to estimate the position on the 
blue side of the Lyman a resonance at which r a b s ~ 1 . Above 
that point Frj — > 1, however this does not represent the main 
behavior of F v for the given assumptions, since also the factor 
1 // v < in a becomes important, so that F v instead actually scales 
like 1 lf v at large x^ (see Sect. 13.4.2] ). 

3.3.3. Effect of line absorption profile on r abs 

It is clear that also the shape of the absorption profile has an 
effect on the frequency dependence of the the absorption op- 
tical depth. As we explained in Sect. 13.21 the effective absorp- 
tion profile, 0* bs is very close the two-photon emission pro- 
file of the 3d-level (see Fig. [2}. For simplicity assuming that 
0* bs = 03d^i s , it is clear that at v > V31 no photons can be ab- 
sorbed in the Lyman a transition, since there 03d^i s = 0. This is 
in stark contrast to the case of a normal Voigt profile, for which 
in principle some photons can be absorbed at arbitrarily large 
frequencies. Considering photons that reach the frequency inter- 
val V21 < v < V31, the fact that there 03d^i s ^ <Pv (see Fig. [2j 
implies that the contribution to the total absorption optical depth 
coming from this region is smaller than in the standard '1 + 1' 
photon formulation. Similarly, at v < V21 the contribution to the 
total absorption optical depth becomes larger than in the stan- 



dard case, because there 03d- 



In Fig. [3] lower panel, we illustrate these effects on r a b s for 
the 3 shell hydrogen atom. However, here we used the full ab- 
sorption profile , 0* bs , which at v > V31 has a small contribution 
from the Voigt profile that is used to model the continuum chan- 
nel (Is *-¥ 2p <-> c). Therefore the optical depths does not vanish 
at v » V3 1 . The additional differences in the values of the optical 
depth seen in Fig. [3] confirm the above statements. Comparing 



with the case for the thermodynamic factor (central panel) it is 
clear that the correction to r a b s due to the shape of the absorption 
profile is not as important. 

One should also mention that setting © a = 1 and f v = 1 we 
obtain the solution F^ as given by Eq. ( 115) . With the comments 
made above, one therefore expects a strong drop in the value of 
Fq for v — > V31, since there r a b s — > 0. Numerically we indeed 
find this behavior (see Sectj4]i. 



3.4. Time and frequency dependence of the effective 
emissivity 

If we look at the definition of a , Eq. ( IScl i. and rewrite it like 



a (z,z') 
(z,z') 



jVem(z')-^' 

N em (z) - Nt 
N em (z) - Nt 



0era(v',z') 



(19a) 
(19b) 
(19c) 



we can clearly see that there are also three sources for the cor- 
rections to the effective emissivity. The first is due to the purely 
time-dependent correction (0 a = ©'), the second due to the ther- 
modynamic correction factor (0 a = 1// V '), and the last due to 
the quantum mechanical asymmet nPl between the emission and 
absorption profile (0 a = 0^). Below we now shortly discuss the 
expected consequences of each of these. 



3.4.1 . Pure time-dependent correction to a 

For a = 0' we are looking at the pure time-dependent cor- 
rection to the emission c oefficient. This correction w as already 
discussed in detail earlier dChluba & Sunvaevl2008bl) . For quasi- 
stationary conditions one would have l = 1. However, in the 
cosmological recombination problem 0' + 1 most of the time. 
This fact leads to significant changes in the shape of the spec- 
tral distortion at different redshifts, where at frequencies x-q < 
only ' + 1 is able to affect the distortion dChluba & Sunvaevl 
l2008bl) . 



3.4.2. Effect of Thermodynamic correction factor in a 



If we only include the correction due to the thermodynamic fac- 
tor f v then we have a = 1 j f v . Since for V — v jfr + V21 one 
has fy + 1, so that due to f v one expects similar effect on the 
shape of the distortion like from l . However, since / v < » 1 at 
large detuning blue-ward of the line center, it turns out that this 
correction can be very large. As mentioned in Sect. 13.31 from 
Eq. (1151 . one naively expects Fq — > 1, but when including the 
factor 1 / fy in a instead we find F v ~ 1 //„ at large xu- 

To show this, it is illustrative to look at the solution for F{v, z) 
in this case, assuming that the quasi-stationary approximation 
(Az/z <k 1 and Av/v <k 1 between the emission and absorption 
redshift of the photons) is possible. Introducing the new variable 

X(v) = £ C bs ( v ') dv '> and usin § T abs ~ T dfe) /v(z) W ~Xl wnere 



10 More clearly here one should refer to the mixture of quantum me- 
chanical processes important for the emission and absorption profile. 
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X' — X( v ')' one has 



f(v,z)=J* S 



*■ ccr r (z')N ls (z') 



H{z'){\+z') 
Jo 



€bs( v '^')e- T * (v ' z ' z) dz' (20a) 

I _ g-Tj fvlXs-Xl 



Mz) 



(20b) 



Since Td/vU's - X\ ~ * 00 f° r x d » 1, there one has F(v,z) — > 
1 / / y (z). In addition one expects F(v, z) ~ 1 // v (z) for xd < 1 . As 
we will show below, for the correction due to the thermodynamic 
factor the scaling F(v,z) ~ l// y (z) indeed is correct at xd » 1 
and xd ~ 1 . However, at xd — > -oo one finds F(v, z) ~ 1 instead. 

3.4.3. Correction due to the quantum-mechanical asymmetry 
between emission and absorption profile 

Since in general the admixture the different transition channels 
to the emission and absorption profile is not identical one does 
expect 0* bs + <p em . For this we can also look at the difference 
between </>* bs and (f> em , which will be given by 



z 



R~ R~ 

2p 2p 



R 



2p 



Uv). 



(21) 



Since in full thermodynamic equilibrium (R 1 ^ ) eq = (^!, p A^2 P ) eq 



2p 



it is clear that {R\~ /fl 2p ) eq = (^p/ ^+ D ) eq , so that Ac/> = 0. On 
the other hand it is known (e.g. see IChluba et alj|2007l) that in 
the cosmological recombination problem R'^~ and R^ + should al- 
ways be very close to their equilibrium values, so that one ex- 
pects <k 1. Only at low redshifts (z < 800) this condition 
may not be fulfilled. However, as we will see below in the con- 
text of CMB power spectrum computations this aspect of the 
problem never becomes significant (see Sect. |4]i. 

4. Changes in the Lyman a spectral distortion 

In this Section we show the detailed dependence of the resulting 
Lyman a spectral distortion on the different corrections that are 
taken into account. As explained above there are three types of 
corrections that are considered here: (i) the time-dependent cor- 
rection, (ii) the thermodynamic correction factor f v , and (iii) the 
dependence on the detailed shape of the effective line emission 
and absorption profile. We start our discussion by first only in- 
cluding these corrections in the computation of r a b s but setting 
& = 1 and 0* = (Sect. 1470. In Sect. @~2] we then also allow 
for 0' + 1, but still set 0^ = 0. Finally, we also include the 
correction due to 0^ + (Sect. 14. 31 ), but this aspect of the prob- 
lem turns out to be not very important. It should be possible to 
deduce all the other combinations from these cases. 

4.1. Corrections related to r abs for & = 1 and & <p = 

As a first case we study the effect of different corrections to the 
absorption optical depth. For this we set 0' = 1 and 0* = 0, 
meaning that in the emission coefficient a we ignore the pure 
time-dependent correction and the one related to the shape of the 
profile. However, depending on the considered case we do allow 
for these corrections in the computation of r a b s . In addition, we 
also discuss the effect of the thermodynamic correction factor, 
f v . but here we include it in both T a b s and a at the same time. 
As explained Sect. 13.4.21 if one would only include f v for r a b S or 
a separately, the corresponding spectral distortion is physically 
not very meaningful. Therefore we omit this case here. 



4.1 .1 . Behavior very close to the line center 

In Fig. [4] as an example we show the Lyman a spectral distortion 
at z = 1100 in the close vicinity of the line center. We com- 
pare the results with the normal quasi-stationary solution (see 
IChluba & Sunvaevl2008bl for details) 



F qs (y,z) = \-e 



(22) 



with x — fo (fy(v')dv'. We show the result obtained for the 
pure time-dependent correctio n to r a b s (dashed curve), w hich 
was already discussed earlier dChluba & Sunvaevll2008bl) . At 

< 4 the distortion is practically identical with the result in 
the quasi-stationary case, while at xd > 4 the time-dependent 
corrections to r a b S start to be important. One can see that there 
F Y = 1 - e~ T » bs > F qs , which as explained in Sect. 13.3.11 is re- 
lated to the fact that including the time-dependent correction, at 
this redshift the value of r a b s is slightly larger than in the quasi- 
stationary approximation, leading tcPI r^ bs > r q bs . 

If we also include the correction due to the shape of the ab- 
sorption profile in the computation of r a b s (dashed-dotted curve), 
then we see that at xd < again the solution is practically iden- 
tical with the solution in the quasi-stationary case. Although one 
does expect some corrections to the exact value of r a b s at dif- 
ferent frequencies below the line centeiP^l since r a b s » 1 the 
effect on the shape of F v ~ 1 - e~ Tl,bs will be exponentially small. 
However, looking at xq > 4 we can see that F v < F qs , imply- 
ing that n < . Since blue-ward of the resonance 0* bs < cpy 
(see right panel of Fig. |2]i it is already expected that the curve 
lies below the one for the pure time-dependent correction (i.e. 
T't 55 t'. )■ But it even turns out that the correction due to the 

abs ~ abs 7 

shape of 0* bs overcompensates the pure time-dependent effect, 
which alone leads to t', > r q ^ . This shows that at the consid- 

abs abs 

ered redshift the correction due to the profile is slightly more 
important than due to time-dependence. 

If we now only include the time-dependent correction and 
the effect of the thermodynamic correction factor (in both r a b s 
and a ) then we obtain the dotted line. As expected from 
Eq. ( |20l , very close to the Doppler core (|xd| < 4) one has 
F v (z) ~ 1 1 'fv(z). We also found this scaling at other redshifts 
(marginally visible in Fig.|5j, as long as the optical depth across 
the Doppler core is much larger that unity, implying that the 
quasi-stationary approximation is valid. However, outside this 
region the distortion differs significantly from the previous cases. 
In particular one finds F v > F qs at xq < -4, which is the result 
of a = 1 lf Y . If we include the thermodynamic correction factor 
only in the computation of r a b s , i.e. setting a = 1, then we find 
F v ~ 1 instead. 

It also turns out that F v < 1 //„ at xd ^ -4. This is in contrast 
to the result F qs (v,z), as given by Eq. < T20b , for which we assume 
quasi-stationary conditions. This implies that in the red damping 
wing deviations from the quasi-stationary assumption become 
important. 

If we in addition include the correction due to the shape of 
the absorption profile, then we can see that at xd < 4 the curve 
basically coincides with the one from the previous case. This 
again is expected since the tiny corrections to the value of a very 
large optical will not lead to visible changes in F v . In addition, 
at xd ^4 one can see that the difference to the previous case 



11 The upper indices will henceforth indicate which correction was 
included, 't' will stand for the time-dependent correction, 'f for the 
thermodynamic factor, and '<p' for the profile correction. When all cor- 
rections are included simultaneously we will use 'a'. 



Since in the red wins i 



(cf. Fig.|2j one expects t* > T q . 
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Fig. 4. Difference in the Lyman a spectral distortion with respect to the quasi-stationary solution in the no redistribution approx- 
imation, F qs (v, z) as given by Eq. (I22K at z = 1100 close to the line center. For all computations shown in the left panel we set 
0' = 1 and 0^ = 0, while in the right we only set 0^ = 0. The cases labeled with / are computed using © a = 1 /fy (left panel) 
and a = ©V/v' (right panel), while for the others we set f v = 1. In addition the quoted correction factors were included in the 
computation of r a b s (v, z s ,z'). We assumed a 3 shell hydrogen atom. 



are about the same as for the difference between the first two 
cases, owing to the fact that the corrections are small and hence 
additive. 

4.1 .2. Behavior at intermediate and large distances from the 
line center 

We now look at the corrections in a slightly wider range around 
the line center. In Fig. [5] we show the same cases as above, but 
now also varying the redshift. As before one can see that differ- 
ences due to the shape of the absorption profile are negligible 
at xu 5 0. Furthermore, on the blue side of the resonance the 
correction due to the shape of the absorption profile is always 
negative, as also seen in the previous paragraph. 

Taking the differences between the curves for a = 1 (first 
two lines) and those for a = 1 // v < (last two lines) one can also 
see that at z = 1200 these are practically the same. However for 
Z = 1000 and z = 800 higher order terms already start to become 
important. For example, at xq = 100 the difference of the curves 
for a = 1 is ~ 2 x 10~ 3 , while it is about 8 x 10~ 4 for those with 
© a = l//v 

If we consider the distortion in an even wider range of fre- 
quencies (Fig.|6]l, then we can make several important observa- 
tion. First, as expected from the discussion related to Eq. < f20b in 
the limit xu — > °o for a = l/fy one finds F v ~ 1 /f v , regardless 
if the normal Voigt profile was used or the effective absorption 
profile, as given by Eq. ( T5bl ). However, for 0* bs = 0v the limit 
Fv ~ l//v is reached closer to the line center than for the ef- 
fective absorption profile. This is expected, since for the 3 shell 
atom the effective absorption profile only has a small admixture 
of the Voigt profile (due to the description of routes connecting 
to the continuum). If <f>*. = (p^ then the limit F v ~ 1 /f v would 
never be reached, simply because at v > V31 the contribution 
to Tabs would be zero. Reducing the admixture of the pure Voigt 
profile therefore moves the transition to F v ~ 1 /f v towards larger 
frequencies. 



The second important observation is that in the frequency 
range V21 < v < V31 on the blue side of the resonance the cor- 
rection due to the shape of the absorption profile is much more 
important than both the pure time-dependent correction to r a b s 
and the correction due to the l// v < factor in a . And finally, in 
the red wing the correction to the spectral distortion is domi- 
nated by the 1 / f v > scaling of a , however, the correction is very 
small, in particular in comparison with the one coming from 0' 
(see Sect.l4~2li. 

4.2. Corrections related to both r abs and & but for 0^ = 

We now want to understand the effect of changes in the ioniza- 
tion history and death probability on the effective emission rate. 
We therefore allow 0' + 1 but still set 0^ = 0. We then again 
discuss different combinations of corrections, like in the previ- 
ous section. As we will see the corrections due to 0' + 1 domi- 
nate at large distances on the red side of the resonance, while the 
shape of the profile is most important for the spectral distortion 
on the blue side of the resonance (see Fig. |6). In the vicinity of 
the resonance basically all the correction factors are important. 

4.2.1 . Behavior very close to the line center 

In Fig.|4]we show the Lyman a spectral distortion atz = 1100 
in the very close vicinity of the line center, now also including 
0' + 1 . If we first look at the curves for f v — 1 (dashed and 
dashed-dotted), then we can see that very close to the resonance 
(\xd\ < 4) the solution is not affected by the inclusion of 0' + 1. 
Due to the huge optical depth across the Doppler core (corre- 
sponding to Az/z ~ 10~ 5 ) there practically any time-dependent 
variation of the effective emission coefficient is erased. However, 
moving towards the wings time-dependent aspects become im- 
portant and in particular now also F v + F qs at xu < 0. 

Also one can see that at the considered redshift the difference 
in comparison to the case 0' = 1 (see Fig. |4]i is very small at 
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Fig. 5. Difference in the Lyman a spectral distortion with respect to the quasi-stationary solution in the no redistribution approxi- 
mation, F qs (v, z), as given by Eq. ( l22l i. at several redshift close to the line center. For the all computations shown in the left column 
we set 0' = 1 and 0* = 0, while in the right we only set 0^ = 0. The cases labeled with / are computed using a = l/f v > (left 
column) and a = 0'//,/ (right column), while for the others we set f Y = 1. In addition the quoted correction factors were included 
for r a bs(v, z$, z0. We assumed a 3 shell hydrogen atom. 



Chluba and Sunyaev: Ly a escape during cosmological hydrogen recombination 



13 



3 shell Hydrogen atom, z = 1200 



3 shell Hydrogen atom, z = 1200 



quasi -stationary 
time dependent T 
time dependent z & $ 
time dependent z & f 
time dependent T & &/ 




'-"lSOOO -10000 -5000 

X D 

3 shell Hydrogen atom, . 



5000 
= 1000 



i|uasi-Mationary 
time dependent T 
time dependent T & (J) 
time dependent T &f 
time dependent T & <j» &/ 




'-^0000 



-10000 10000 

3 shell Hydrogen atom, % = 800 



20000 



quasi -stationary 
time dependent z 
time dependent z & $ 
time dependent T & / 
time dependent T & <j» &/ 




-25000 



I I I I I I I L 

-15000 -5000 5000 




! ?15000 



quasi-stationary 
time dependent T & 0, 
time dependent T & & (ft 
time dependent T & 0, & / 
time dependent T & 0, & $ &/ 

i I i , I I 



3 shell Hydrogen atom, ; = 1000 



10 

10° 
10"' 




Maximum of the 
Ly Ct distortion 



'Soooo 



quasi-stationary 
time dependent T & 0, 
time dependent T & t & *t> 
time dependent T & t & / 
time dependent T & 0, & (i) &/ 

i I I , i I 



-10000 



3 shell Hydrogen atom. 




20000 



10 

10 1 
10° 
L0- ! 



15000 



25(1011 



10' 

io-' 

10"' 

io' 1 

10"' 
10"" 




quasi-stationary 
time dependent z & t 
time dependent T & © t & <fi 
time dependent T & © t & f 
time dependent T & © t & <J> & f 



Fig. 6. The Lyman a spectral distortion at different redshifts and in a wide range of frequencies around the line center. For the all 
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-*d ^ -4. There clearly the correction due to the thermodynamic 
factor f v (dotted and solid lines) is dominant. However, at xu > 4 
the time-dependent changes in the effective emission coefficient 
lead to a correction that is similarly important as the one due to 
f v . One can also see that all the corrections add roughly linearly. 
Note that at xq < 4 the curves are not affected when accounting 
for the corrections to the shape of the line profile. 



4.2.2. Behavior at intermediate and large distances from the 
line center 

Looking at the right column of Fig. [5]we can see that at xq < 
the correction due to the inclusion of /„ dominates at high red- 
shifts, while at z ~ 1000 the correction due to 0' is already com- 
parable, and clearly dominates at low redshifts. In addition, at 
xd Z one can see that at high redshift all the corrections practi- 
cally add linearly, while for z = 1000 and z = 800 the correction 
due to the inclusion of l practically disappears when including 
the correct shape of the effective absorption profile. Also when 
including the thermodynamic correction factor the large excess 
of photons seen for the case f v — 1 and a b s = <pv practically 
vanishes. This implies that the sel f-feedback effect at low re d- 
shifts that was reported elsewhere (IChluba & Sun yaev 20071) is 
expected to disappear. As explained there, this unphysical aspect 
of the solution in the '1 + 1' photon pictures is due to the fact that 
the Voigt profile in principle extends up to arbitrarily large fre- 
quencies, so that photon emitted at z ~ 1400 in the very distant 
blue damping wing will still be able to reach the line center at 
low redshift, strongly enhancing the photon occupation number. 
However, when including the thermodynamic correction factor, 
due to the exponential enhancement of the absorption optical 
depth at large distances blue-ward of the resonance, such pho- 
tons disappear much before this. Similarly, when including the 
shape of the effective emission profile such excess of photons 
will never be produced in the first place, so that from this region 
the residual correction due to 0' + 1 is much smaller. 

Looking at the spectral distortion in a very large range of 
frequencies around the line center (Fig. [6} it is clear that at all 
redshift the shape of the distortion is dominated by the correction 
due to 0' + 1 for xu < 0. Both the thermodynamic factor and the 
shape of the absorption and emission profile only lead to small 
additional modifications there. The largest correction is due to 
the fact that the 3d two-photon emission profile is larger than 
the Voigt profile at v — > vn/2, explaining the small addition of 
photons in comparison to the case 0* bs = (py seen close to the 
lowest frequencies shown in the figures. On the other hand, at 
frequencies above the line center clearly the correction due to 
the shape of the line profile is most important. In the line center 
all sources of correction contribute to changes of the Lyman a 
spectral distortion with respect to the quasi-stationary solution. 

4.3. Corrections related to r a b s and 0' including 0* + 

We also ran cases for 0** + 0. However, the correction was al- 
ways extremely small. Therefore we decided to omit the corre- 
sponding plots for F v . As mentioned in Sect. 13. 21 this is expected 
since the deviations of R'^/R^ and R'^/R^p from their equilib- 
rium values is always very small in the relevant redshift range, so 
that R'£ /R+ p * tfjp/^Ip' and hence Cbs ~ 0em- In Fig. Q we ex- 
plicitly show this fact. As an example, for the 3d-channel one can 
see that at z ~ 1 100 one has Rl d '~/R; -Rl d - + /RL ~ 10 6 - 10 5 . 



10 shell Hydrogen atom 



Therefore one would expect corrections to the effective escape 
probability at the level of ~ 10~ 4 - 10~ 3 %, which is clearly neg- 
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Fig. 7. Source of the asymmetry between the absorption and 
emission profiles. We present A(z) = R'^/^p ~ ^'ip^ip ^ or a 
given level as a function of redshift. For all curves we used the 
solution for the 10 shell hydrogen atom. 



ligible in the context of the CMB anisotropics. We confirmed 
this statement numerically. 

One should also mention that although at low redshifts the 
expected difference between emission and absorption profile in- 
creases, there the value of the effective escape probability is 
dominated by the contribution from the Doppler core, where de- 
tails of the profiles will not matter very much. In addition correc- 
tions to the escape probability will not propagate very much into 
the ionization history, so that here we do not discuss this point 
any further. 

5. Changes in the effective escape probability 

5.1. Effect of the thermodynamic correction factor 

First we consider the normal Lyman a transfer equation ([TJ, but 
including the thermodynamic correction factor f v , in order to 
correct for the small imbalance in the emission and absorption 
process in the line wings coming from the normal '1 + 1' photon 
formulation (see Sect. 12.1.11 ). In this case one has 0;(v) = <pv(v) 
and hence em = </>* bs = 4>v(v), so that from d8bl one can find 



F f Y = J & f ( Z , z') d z ,e- T -^ v ' z '- z) dz' 



& f (z,z) = 



(23a) 
(23b) 



with T a b s as given by Eq. d5al l but for 0* bs (v) — > <f>v(v)- 

If here one drops the factors due to f Y in the definition of 
f and also r a b s , one obtains the purely time-dependent case , 
F\,, that has been addressed earlier dChluba & Sunvaevll2008bl) . 
However, here the term 0' naturally is given by the line center 
values of N em (z') - Nf, which in the previous work had to be 
enforced for consistency with the standard approxim ations (see 



enrorcea ror consistency witn tne standard approximal 
comments in Sect. 3.4.1 of lChluba & Sunvaevll2008bl) . 
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Fig. 8. The functions Tl = f v F f v at redshift z = 1000. In the left panel we show the difference of AA = Tl - Ff, while in the right 
we directly show A = T v - The cases labeled with / are computed using a = & I f r , the others with a = l . In addition the quoted 
correction factors were included in the computation of T a b s (v, z s ,z')- For comparison in the right panel we also directly show the 
normal quasi-stationary solution f\ qs (double dash dotted line) and Fy (dotted line) that includes the correction due to the profile 
0* bs + 0v (see Sect. 15.2.2b . We assumed a 3 shell hydrogen atom. 



To understand the role of f v in the final correction to P e g, 
Eq. ( |12bt , we first look at the term APj n( j, Eq. (|12ct . It is clear 
that for <f> em = 0* bs = <py one has APi nc j = 0, since in this case 
Aijo em = A^* bs = 0. As we will show below (Sect. 15.2. IK AP; n d + 
when taking into account the effect of stimulated emission in the 
blackbody radiation field, however, the contribution to the final 
correction turns out to be negligible. 

If we now look at the definition of P, Eq. ( llObt , then for the 
considered case we can rewrite it as 

P s = j <p y (v)[l-Tl]dv (24) 

where we introduced the function 

7?=Mz)F f Y . (25) 

This representation allows to directly reveal the expected differ- 
ences in comparison with the standard quasi-stationary case, for 
which one has (cp em = i^* bs = ipy and f v = 1) 

TT =Ff = l- e- Ti e TiX . (26) 

If instead of Tl we insert this into Eq. ( f24l , we directly obtain 
Fd = 1 ~T d Td • Therefore we can write the correction with respect 
to the quasi-stationary solution as 

API = P f - Pi = - J Mv) [Tl - T v qs ] dv. (27) 

This expression shows that for the correction to the effective es- 
cape probability it in fact is not important how F f v behaves, but 
how Tl looks. Since f v is a very strong function of frequency 
this makes a big difference in the conclusions, as we explain be- 
low. Furthermore, any difference to the value of TT will lead 
to a contribution to AF^ with opposite sign. Also, it is possible 
to compute the partial contributions to the total correction in the 
escape probability by only integrating AP f d as defined by Eq. ( f2Tb 
over a given range of frequencies. 



5.1 .1 . Behavior of Tl and the expected corrections to F eff 

In Fig. [8] as an example we show Tl at redshift z = 1000 (dotted 
lines). For comparison we also show the pure time-dependent so- 
lution (dashed curves) for which Tl = F' v . Note that we include 
the correction terms in both r a b s and a . 

If we focus on the behavior at -100 < a'd < -4, then we 
can see that although in the considered cases F\ > F\ > Ff 
(compare with Fig. [5] right column, case z = 1000), for Tl 
one has Tl > T^ > Tl- This change in the sequence is only 
due to the factor /„ < 1 in the definition of Tl, which appears 
due to a b s = fv</>Zh S m th e esca P e integral, Eq. dlObl i. Instead 
of an additional negative contribution to the escape probabil- 
ity (AP f < AP' < 0), as it would follow from the differences 
seen in F v , when including the thermodynamic correction factor 
f v one therefore expects a positive contribution from the con- 
sidered spectral region. This effect becomes even more appar- 
ent when looking at the very distant red wing, where clearly 
Tl > TT » Ty, owing to the exponential cutoff introduced 
by fv The behavior shows that in the very distant red wing the 
excess Lyman a photons no longer supports the flow of electrons 
toward higher levels and the continuum. These photons only un- 
dergo line scattering events, with tiny shifts in the frequency due 
to the Doppler motion of the atom, but do not die anymore, and 
therefore have effectively already escaped. As we will demon- 
strate below (Sect. 15.1.41 and |5.1.5t in fact the main correction 
due to f v is coming from the change in the absorption cross sec- 
tion in the red wing of the line profile. The real modifications of 
the spectral distortion due to f v are not so important. 

If we now look at the behavior in the range 4 < < 100 
we can see that the sequence F v > F f v > Fy^ does not change 
when considering Tl > Tl > Ty S , but Tl becomes more sim- 
ilar to Tl- From the behavior at -100 < < -4 and the 
strength of the changes seen there one might have expected that 
at 4 < x D < 100 also Tl > Tl, since |A/ y // v | in fact is sim- 
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Fig. 9. Relative difference in the effective escape probability with 
respect to the Sobolev escape probability: effect of the thermo- 
dynamic correction factor. 



ilar in both regions. However, at 4 < xu < 100 the spectral 
distortion F v is very steep, so that small changes A/ v // v can- 
not affect the shape of T$ so much in comparison to F f v . Only 
when going to much larger distances on the blue side of the 
resonance, where f v again behaves exponentially, one can see 
Ty » Ty > TT although there F' v » F f v (compare with Fig.|6] 
right column, case z = 1000). In comparison to the pure time- 
dependent correction from the range 4 < < 100 one therefore 
expects a slightly smaller (negative) correction to the total value 
of P, while the contributions from very large distances in the 
blue wing should be significantly larger than in the pure time- 
dependent case. However, here it is important to mention that 
these very distant wing contributions will always be very minor, 
since the Voigt profile drops like oc 1 /x^ (see Sect. 15.1.31 ). 

Finally, looking at the central region -4 < jcq < 4 we can see 
that there now « Ty ~ '. This is because F s v « l/f v in that 
regions (see Sect. 14.1. It , so that Tl ~ 1 with very high accuracy. 
This implies that although the Lyman a spectral distortions in the 
case F\ and F"y S look rather different the correction factor f v does 
not lead to any real correction in the escape probability from in- 
side the Doppler core. There everything is well described within 
the quasi-stationary assumption, for which the whole Doppler 
core reaches thermodynamic equilibrium with the ambient radi- 
ation field, but now also including the small additional variation 
of the photon distribution over v. 

5.1.2. Correction to the escape probability 

In Fig. [9] we show the result for the effective escape probability 
and the effect of f v . For comparison we also show the result for 
the pure time-dependent co rrection (dotted line) which was al- 
ready discussed elsewhere (IChluba & Sunvaevll2008bl) . At low 
redshifts we indicate the rise in the amplitude of the correction, 
which was attributed to the late self feedback of Lyman a pho- 
tons for this case. For the other cases we also point towards the 
difference which arises due to the no scattering approximation. 
It is due to the differences in and Ps itself which close to 
the maximum of the Thomson visibility function at z ~ 1100 



are negligible, but become notable both at very low and very 
high redshifts ( IChluba & Sunvaevll2008bl) . However, there the 
changes have no effect on the free electron fraction. 

When now including the thermodynamic correction factor in 
the computation of the absorption optical depth and the effective 
emission rate a , but setting 0' = 1 (dashed line), one can see 
that the correction to the escape probability becomes positive at 
all redshifts, with a maximum of AP/P ~ 11% at z ~ 1350. As 
we will see in Sect. 15 . 1 .41 and I5TT31 bulk of this total correction 
is coming from changes in the absorption process on the red side 
of the resonance, where in this case F v =s Fy S » 1 (cf. Fig.|5]and 
[6] left column). 

When also including the variation of 0' + 1 over time (solid 
line) the result changes significantly, shifting the maximum of 
the correction AP/P ~ 10% to z ~ 1450, but still leading to 
AP/P > everywhere. However, especially the low redshiftpart 
is strongly modified, reducing the total correction by a factor ~ 2 
at z ~ 1 100. We also show the infered correction due to l alone, 
which was obtained by taking the difference between the curves 
labeled (I) and (II). The result shows that the final correction 
close to the maximum of the Thomson visibility function has im- 
portant contributions from both f v + 1 and the time-dependence 
of the problem. 

5.1 .3. Effect at large distances blue-ward of the line center 
and the Lyman a self feedback 

The thermodynamic factor clearly strongly changes the behavior 
of the correction to the effective escape probability. The purely 
time-dependent correction is no longer dominant, and in partic- 
ular the thermodynamic factor removes the self feedback prob- 
lem of Lyman a at low redshifts, which was already realized 
to be an artifact of the st andard '1 + 1' photon formulation 
(IChluba & Sunvaevll2008bl) . This is due to the fact that when 
taking f v into account, exponentially fewer photons remain in 
the photon distribution at large distances on the blue side of the 
resonance (cf. Fig.|6]where F v ~ 1 // v at large jcd)- As explained 
in Sect. [3] every photon emitted at xo » 1 will be re-absorbed 
quasi-instantaneously. This is due to the exponentially larger 
amount of CMB photons red-ward of the the Balmer a line, than 
close to the line center, so that line absorption is more effective. 
The main process for the death here is the Is — > 2p — > 3d, where 
the last step is considered to lead to a complete redistribution, 
so that the absorbed Lyman a photon (Is — > 2p) will most likely 
reappear close to the line center. Note that in this Section we still 
have not included the two-photon corrections to the shape of the 
absorption profile, but already now the two-photon character of 
the process leads to this conclusion. 

However, one has to mention that basically all the photons 
present at these large distances contribute to the escape integral. 
This is b ecause Ty ~ 1 while TT <sc 1, so that J* - TT ~ 1 
(see Sect. 15.1. H and Fig. [SJ. Therefore in this case the very distant 
blue wing contribution to the value of AP\ behaves like 



API 

DW 

~p7 



4>m dv' 



7T X \ 



bw ' 



(28) 



at z 



~ 1100 implying an additional < 0.1% correction from 
> 10 4 , showing that the absolute contribution becomes 
negligible far beyond that point. In numerical computations we 
therefore typically integrated the profile up to xu ~ 1.5 x 10 4 , in- 
troducing a very small error for cases with f v + 1 and practically 
no error when also including the correct shape of the emission 
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Fig. 10. Relative difference in the effective escape probability 
with respect to the Sobolev escape probability: effect of the ther- 
modynamic correction factor in the blue and red wing of the 
resonance. For the upper panel we set 0' = 1 and 0* = 0. The 
lower panel was computed also including 0' + 1 . 



profile, since then naturally very few photons appear at larger 
distances blue-ward of the resonance (e.g. see Fig. [8] solid line). 



5.1 .4. Contributions from the blue and red wing 

Returning to the correction at higher redshifts (z > 800 - 900), it 
was already sho wn earlier, that there the Lym an a self feedback 
is not important (IChluba & Sunvaevll2008bl) . From Fig. |9] it is 
clear that both considered cases for the effect of f v imply that 
at a given redshift effectively fewer photons support the flow of 
electrons towards higher levels and the continuum than in the 
quasi-stationary case, albeit the fact the more photons are pro- 
duced. The latter statement can also be confirmed by looking at 
Fig. [6] and the amplitude of the Lyman a distortion in the dis- 



tant red wing around its maximum. Note that these photons have 
already been emitted at z ~ 1400. 

However, where does the main correction to the escape prob- 
ability come from at these redshifts? Looking at Fig. [8] we can 
see that at < xd < 100 the function Ty is very similar to 
Tl, which results in the pure time-dependent correction. Also at 
slightly larger distances (100 < xd 55 10 3 ) one still has Tl ~ Tl- 
Therefore, one does not expect very dramatic changes in the con- 
tribution to the effective escape from this part of the Lyman a 
distortion in comparison to the purely time-dependent case. 

On the other hand, in the red wing one finds Tl < Tl, and 
at very large distances one even has Tl <K 7^. Physically this 
reflects the fact that due to f v the re-absorption process in the 
distant red wing is exponentially suppressed, so that there pho- 
tons can escape more directly, than in the normal '1 + 1' pho- 
ton formulation. This is now related to the exponentially smaller 
amount of CMB photons blue-ward of the Balmer a line, so that 
the main absorption channel Is — > 2p — > 3d becomes practically 
inactive at xu < -10 3 . It therefore is the modifications in the 
red wing absorption process from which one expects the largest 
effect in connection with f v at high redshifts, before the appear- 
ance of the self feedback problem. 

Numerically we studied this statement by simply assuming 
that at < the solution is given by the quasi-stationary result 
(implying Tl = T^ ~ 1), while at < xd we used the real so- 
lution for T v in the considered case. In this way it is possible to 
separate the 'blue wing' contribution to the total correction in the 
effective escape probability, and similarly one can obtain the 'red 
wing' contribution. In Fig.[10]we show the comparison of these 
computations for the cases 0' = 1 and 0^ = (upper panel; only 
the corrections to t are included and a = l// v <), and 0' + 1 
and 0* = (lower panel; a = & lf v >). For comparison we also 
show the results obtained for th e purely time-dependent co rrec- 
tion in the considered cases (cf. IChluba & Sunvaevll2008bl) . 

In the first case (upper panel), one can clearly see that the 
blue wing contribution from Tl is very close to the purely time- 
dependent result (dotted curve), which itself in fact has no sig- 
nificant contribution from the red wing in the first place (dia- 
monds). This shows that for this case the effect of f v is not im- 
portant at xd > 0. One can see that indeed the main correction 
due to the effect of f v arises from the red wing, and that this 
correction is significantly larger than the time-dependent case 
alone. 

If we look at the comparison in the full time-dependent 
case (lower panel), one can see that when including the cor- 
rection factor f v , at high redshifts the blue wing contribution 
(dashed curve) is a bout 50 % of the total result presented in 
IChluba & Sunvaevl d2008bl) . At high redshifts the blue-wing 
contribution in the pure time-dependent case (stars) practically 
coincides with the one that includes f v , implying again that the 
blue wing contribution is not affected much by the thermody- 
namic correction factor. However, one can see that at low red- 
shifts f v is very important for avoiding the self feedback prob- 
lem, as explained in Sect. 15.1.31 Note tha t in contrast to the 
curve q uoted 'AP X (only blue)' in Fig. 8 of IChluba & Sunvaevl 
(I2008M) . here the blue wing contribution takes into account the 
time-dependent correction on t and 0' simultaneously. 

Looking at the red wing contributions for this case one can 
see that for f Y = 1 (diamonds) now t he contribution is non zero . 
This was also already seen earlier (IChluba & Sunvaevl l2008bb 
and is due to the fact that Fj, + 1. However, the contribution 
from the red wing is much larger when including f v + 1, and 
in particular it is only positive due to the fact that Tl < 1 at all 
frequencies, so that AP > 0. The conclusion clearly is that the 
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Fig. 11. Relative difference in the effective escape probability 
with respect to the Sobolev escape probability: effect of the ther- 
modynamic correction factor at different distances to the line 
center. For the upper panel we set 0' = 1 and = 0. The 
lower panel was computed also including ©' + 1 . 



dominant correction due to the inclusion of f v is coming from 
the red wing of the Lyman a resonance. 

5.1 .5. Simple estimate for the red wing correction 

We can also perform another rough estimate for the expected 
correction, assuming that in the red wing Tl ~ f v as suggested 
by Fig. [8] This will overestimate the result, since with the inclu- 
sion of f Y alone one already obtains Tl > f v . In comparison with 
the quasi-stationary case (TT ~ 1) we then have 
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1100 



(30) 



At redshift z ~ 1 100 one has ~ 6.8 x 10 4 , so that only from the 



dx D . (29) 



AP 

red wing correction one expects ~ 11%, while at redshift 

'& 

z ~ 1350 one finds ~ 16%. Looking at the upper panel in 
Fig- [TOjshows that this is the right order of magnitude, although 
the final correction is about 1 .5 times smaller than given by this 
simple estimate. 

5.1 .6. Dependence on the distance to the line center 

Finally we want to look at the dependence of the correction on 
the distance to the line center. For this we computed the results 
including the deviation for the quasi-stationary case in a given 
range around the line center. The results of these computations 
are shown in Fig. QT| Clearly a large fraction of the total correc- 
tion is coming from large distances (10 < |xdI 5 100 - 1000) 
from the line center, while the contributions from within the 
Doppler core (|xd| < 4) are very small. The latter result again 
reflects the fact that there neither f v deviates very strongly from 
unity, nor does any time-dependent effect become important. 
The Doppler co re can be considered quasi -stationary (for more 
explanation see IChluba & Sun yaev 2008b) and in full equilib- 
rium with the line center value. 

5.2. Dependence on the shape of the absorption profile 

As a next step we want to understand how the two-photon cor- 
rections to the shape of the effective line profile affect the escape 
probability. For this we ran computations only including the fact 
that 0* bs + (f>\, but for the moment we shall neglect the correc- 
tion due to f v and also assume 0* bs = <f> em inside a . As explained 
in Sect. |4.3| the latter correction for our purpose is negligible dur- 
ing cosmological recombination, but the inclusion of f v + 1 is 
still expected to be very important, as we will discuss in Sect. 15.31 
for the combined effect. 

5.2.1 . Correction due to AP ind 

Looking at Eq. dl2ct , it is clear that the contribution AfVd to 
P e ff is purely due to induced effects, since in vacuum one would 
find Aip em = Ai£* bs = 0, and hence AP; n d = 0. In Fig. Q~2] we 
present the deviation in the normalization of the emission and 
absorption profile from unity, which have been computed using 
Eq. (|6]l and d5bl ). First one can see that at basically all redshifts of 
interest in the recombination problem Aip em m A(jo* bs , implying 
that APj n d ~ Ai£ em . Since A(p em < 6 x 10~ 9 at all shown redshifts, 
comparing AP; IK j with Pd implies that the associated correction 
APs should never exceed the level of 0.1%. In a more detailed 
computation we find a maximal correction of APs/Ps ~ 0.06% 
at z ~ 1300. In addition, this correction is practically canceled 
by another contribution, APi n d, n0 rm = _ A^* bs , which appears as a 
result of stimulated emission on the overall normalization of the 
effective absorption profile (see Sect. 15.2.21 ). We therefore can 
neglect this term in the following. 

We would like to mention that the main contribution to 
<Pem ^ 1 is coming from the region v ~ V31. This can be seen 
in Fig. 12 where only there the effective two-photon emission 
profile differs significantly from the vacuum profile, ^* mvac . At 
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Fig. 12. Deviation of the normalization of the different line- 
profiles from unity. The effect of stimulated emission in the am- 
bient CMB blackbody radiation field was included. 



this large distance from the line center the overall profile already 
dropped by a factor of 10 11 - 10 12 relative to the line center (the 
value there is ~ 1/ yjn ~ 0.56). From Fig. [2] one can see that 
there y>* m - <p* em vac ~ 10" 12 over a region Axq ~ 10 3 , so that one 
expects APi n d ~ 10~ 9 at z ~ 1300, which is in good agreement 
with the full numerical result. 



5.2.2. Expression for P eff and its connection to F v . 

In order to understand the corrections in the effective escape 
probability and its relation to the differences in the spectral dis- 
tortion we again look at the definition of P, Eq. ( llObb . with 
f v = 1 and rewrite it as 

P t4 = f[<Pv - Vlj dv + J ^ [1 - F l f] dv. (31) 
The first integral is given by 

AP ind ,norm = J [<fv ~ ^s.vaJ dv + J fes.vao ~ Vabs] dv 



abs' 



(32) 



Note that f[<pv -y* bs vac ] dv = even though the partial contribu- 
tions from the red and blue wing are non zero. As mentioned in 
Sect. 15.2.11 APind^orm cancels with the correction due to AP; n d, 
so that we finally have 
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.[l-^dv 



(33) 



This expression now allows to compute the correction to the es- 
cape probability. 

In order to find out how the shape of the profile enters into 
the problem, it is illustrative to look at the result for when as- 
suming quasi-stationary conditions, but including the correction 
due to the profile. In this case one has t^ s = Td(z) <^* bs dv' 
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Fig. 13. Relative difference in the effective escape probability 
with respect to the Sobolev escape probabili ty: effect of the line 
profile . For the curves labeled with 'CS09' dChluba & Sunvaevl 
l2008bl) we used the standard Voigt profile, while for the others 
we included the two-photon corrections for the 3s and 3d chan- 
nels. 



and with Fy' qs (v) = 1 - e resulting from Eq. ( fTBI l one finds 



a x ; * p^ (34) 



Here we have introduced the variable xl( v ) — f Vabs dv ' and 
Xlco - Xl(°°)- This expression shows that in the quasi-stationary 
case the shape of the pr° n l e does not matter. The result will still 
be extremely close tq3 Pd- This also implies that the shape of 
the profile is only entering as correction to correction, i.e. com- 
bination of time-dependence and profile. Therefore the changes 
due to the exact shape of the profile are expected to be smaller 
than the corrections due to f v and the time-dependence alone. As 
we will see below they still are not negligible, in particular when 
taking all corrections into account simultaneously (Sect. |5~3l . 

For the correction to the effective escape probability we can 
therefore finally write 



AF* 



P ',<t> 

eff 



Pd 



Idv. 



(35) 



Now it is clear that for the contributions of the total correction it 
is important how F*f deviates from Fy' qs . For F^f < Ff' qs one 
will have a positive contribution and for Fy > ,Fy' qs a negative. 
In addition due to the appearance of 0* bs in the outer integral, the 
red wing contribution (0* bs > <py) will be slightly overweighted, 
while the blue wing (<p* dbs 55 0v) will be underweighted. 

We would like to mention that Ff' qs (v) behaves very similar 
to P qs (v). In particular at frequencies xq < 4 it also becomes 
very close to unity, so that there F"*' qs (v) * F qh (v). The main 
difference appears in the blue wing of the line, where pf' qs (v) 
depends strongly on the differences between 0* bs and <py. Both 
aspects can be seen in the right panel of Fig. [8] 



13 There is a tiny difference due to the fact that^* , 
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Fig. 14. Relative difference in the effective escape probability 
with respect to the Sobolev escape probability: correction due to 
the shape of the absorption profile in the blue and red wing of the 
resonance. For the upper panel we set 0' = 1 and ®* = 0. The 
lower panel was computed also including ©' + 1 . In all cases 
we used f v = 1. Those curves labeled with <p were computed 
including the 3s and 3d two-photon corrections. 



5.2.3. Total corrections and the blue and red wing 
contributions 

In Fig.[13]we show the results for the total correction to the effec- 
tive escape probability when including the modifications in the 
shape of the absorption profile. We used the expression Eq. d33l 
to compute the different curves. Again for comparison we also 
give the results for the time-dependent corrections only. As one 
can see the main effect of the profile is the removal of the self 
feedback at low redshift. Other than that in the considered cases 
the modifications in comparison with the time-dependent result 
are rather small (less than ~ 10% of AP/P for the cases with 
0< = 1 and less than ~ 1% of AP/P for those with © l + 1). 



This shows that there is a cancelation of corrections from the 
red and the blue wing, since it is clear that already the modula- 
tion of these partial contributions due to the presence of 0* bs in 
the outer integral of Eq. ( 1331 should have some effect, even if it 
is of higher order in the correction. To understand the results we 
therefore look at the differential contribution from the red and 
blue wing separately. For Fig. [14] we ran computations includ- 
ing the corrections to the quasi-stationary result either on the red 
or blue side of the resonance. We compare the total and partial 
corrections in both the pure time-dependent case and when also 
including the shape of the absorption profile. 

Corrections in the case of ©' = 1 

When only including the profile corrections to r a b s (upper panel, 
©' = 1 and ©^ = 0), we can see that the effect of </>* bs is not lead- 
ing to any important correction from the red wing. Looking at 
Fig.|5]and|6](left column), it is clear that in there Fy « Fy S « 1 . 
Since also F" qs * Fy ' qs (see Sect. I5.2.21 i one has Fy - Fy w 
and hence with Eq. (l33T > one expects a partial contribution of 
AP*j ~ 0, confirming the above behavior. The small positive 
bump seen at z ~ 1400 is mainly due to the fact that we started 
our computation of the spectral distortion at some particular time 
in the past (z s = 2000), so that in the considered case the nu- 
merical solution for the spectral distortion, as computed using 
Eq. (l8bl i. drops towards zero below some distance xu < -10 4 in- 
stead of staying close to unity. Therefore we have Fy - ff' qs « 
-1 and hence AP^(x D < -10 4 ) ~ j 0* bs dx D ~ few x 10~ 3 from 
that region, explaining this small excess with respect to P d - We 
also checked this statement numerically by increasing z s as ex- 
pected finding that the bump became smaller. When also includ- 
ing the time-dependence of the emission coefficient (©' ^ 1) this 
small inconsistency of our computation is no longer important, 
since the spectral distortion by itself drops very fast toward zero 
(cf. Fig. [5] and [6] right column). 

Again looking at the upper panel in Fig. [14] we can also 
see that the largest contribution to the total correction is coming 
from the blue wing, and that the difference to the time-dependent 
case is rather small, with AP^f being slightly smaller. This can 
be understood when writing rf = rVf + At 1 "* with 

Ar d (z')^ bs (v',z')dv' 

= 1 AT d (z')^v(v',z)dv' + I AT d (z')A^* bs (v',z')dv'. 

(36) 

with the abbreviations 1 + z' = ^(1 +z), Ard(z') = Td(z') - r<i(z), 
and A^* bs (v',z') = y?* bs O' ', z') - <py{v',z). This approximation 
shows that one has At 1 '* « f Ata(z!) <fy(v' ,z) Av' = At', since 

the profile correction with respect to r^ s is already of higher 
order. For a = 1 one therefore expects 



AP 
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AT ' -l]dv 
At ' -l]dv = AP< 



(37) 



confirming the result seen in Fig.[14]in lowest order. 

The remaining difference is mainly due to the second order 
term in Eq. £36} At^ - At 1 = J^°AT d (z') A^* bs (v',z')dv', which 

we neglected in Eq. (f3Tb . The factor ip* hs e~ T ^ plays a minor 
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role, since the function e _AT <s - 1 is varying much faster. We 
confirmed these statements numerically, finding that for the con- 
sidered case the modulation of the blue wing correction due to 
0* bs in the outer integral can be neglected. 

Corrections in the case of ©' + 1 

To understand the result in the case when also including the 
change in the effective emission coefficient 0' + 1, we again 
look at the red and blue wing contribution separately. Since in 
the red wing Fy ' qs * 1 and because we already saw in Sect. 14.21 
that there the solution for the spectral distortion is dominated by 
the correction due to ©' + 1 only (cf. Fig.|5]and|6]right column), 
we expect that the partial contribution form the red wing will be 
very close to 



3 shell Hydrogen atom 



dv, 



(38) 



where F\ is the spectral distortion for the purely time-dependent 
case. Therefore to lowest order again one will have APf w * 
AP[. W , but in next order this correction will in addition be slightly 
larger in amplitude due to the fact that for xq < one has 
<p* ihs > (py. We can see that this indeed is true comparing the dia- 
monds with the dash-dotted curve in the lower panel of Fig. [14] 

For the blue wing one can argue in a very similar way as 
above. We know that for & = 1 the correction to the escape 
probability is basically given by the time-dependent correction 
in the value of r a bs, but only to higher order due to </>* bs + <py. 
If now including t + 1 the lowest order correction will still be 
given by the purely time-dependent case. The additional mod- 
ulation of the resulting spectral distortion by 0* bs < <py will in 
addition lead to a small decrease in the total amplitude of the 
contribution to the correction. Again this can be seen in Fig. FBI 
when comparing the stars with the short-dashed curve. Only at 
low redshifts the shape of the profile determines the amplitude 
of the correction, removing the self-feedback problem. This is 
because unlike in the case of 0v photons are emitted only in a 
limited range of frequencies. This avoids that photons which are 
released at z ~ 1400 and xd » 10 4 will redshift into the Lyman 
a line at z < 1000, as seen in the norma l ' 1 + 1 ' formulation of 
the problem dChluba & Sunvaevll2008bl) . 

Furthermore, it is clear that the sum of both the red and blue 
wing contribution should again be very close to the purely time- 
dependent case, since the modulation of the contributions form 
the red (enhancement) and blue wing (suppression) in lowest or- 
der will cancel, due to the symmetry around the line center. 



5.3. Combined effect of f v and < 



abs 



With the derivations in the previous Sections it is now straight- 
forward to understand the results for the combined effect of all 
corrections. Following the same line of thoughts we obtain 



[f v Fl-Ff' qs ]dv. 



(39) 



where in F* we include all the corrections simultaneously. 

As in the previous Section it is now clear that on the red side 
of the resonance the profile will enter the computation mainly 
due to its presence in the outer integral. For Ff, and xd < it only 
leads to a very small correction (cf. Fig. [5] and [6] right column). 
However, on the blue side of the resonance the profile correction 
again can be neglected in the outer integral, but should be taken 
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Fig. 15. Relative difference in the effective escape probability 
with respect to the Sobolev escape probability: combined effect 
of the thermodynamic correction factor and the shape of the line 
profile. 



into account when computing the difference A v = f v F^ - ff' q \ 
Also one can conclude that the shape of the profile plays the key 
role in removing the low redshift self-feedback problem. The 
latter statement can be confirmed when looking at the shape of 
the Lyman a distortion at intermediate to high frequencies blue- 
ward of the resonance (cf. Fig. [5] and [6] right column), which is 
clearly dominated by the profile rather than /„. 

Therefore in lowest order one expects the total correction to 
be the superposition of the time-dependent correction and the 
one from the thermodynamic correction factor, where on the red 
side of the resonance each of them is modulated by the shape 
of the profile in the outer integral of Eq. ( |39l in addition, while 
on the blue side the contribution is slightly suppressed due the 
profile corrections to A„. Here it is important that because f v 
strongly changes the symmetry of the problem (the main correc- 
tion is coming from the red wing as shown in Sect. 15. Il l, it is 
clear that the main effect of 0* bs + (py will be an enhancement of 
the final correction. 

In Fig.[l5]we present the results from our numerical calcula- 
tion for different cases. Indeed we find that when including the 
shape of the profile the corresponding correction is slightly in- 
creased as explained above. 



5.3.1 . Dependence on the included number of shells 

For the purely time-dependent correction it has been shown that 
in particular at low redshifts the result depends strongly on the 
total number of shells that were included in the computation. 
Since there the correction was very strongly dominated by the 
self-feedback of Lyman a photons, here we do expect this de- 
pendence on the number of shells to be more mild. 

In Fig. [16] we show the results of our computations for 3, 5 
and 10 shells. The changes between the 3 and 5 shell cases is still 
rather significant, but the difference between the 5 and 10 shell 
case is already very minor. This show that in our description the 
total correction is already converged when including ~ 5 shells 
into the computation. 
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Fig. 16. Total relative difference in the effective escape probabil- 
ity with respect to the Sobolev escape probability: dependence 
on the total number of shells. Note that the curves for the 5 shell 
and 10 shell cases practically coincide. 
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Fig. 18. Changes in the free electron fraction: separate contribu- 
tions due to the time-dependent correction, the thermodynamic 
factor and the shape of the profile. 
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Fig. 17. Relative difference in the effective escape probability 
with respect to the Sobolev escape probability: separate con- 
tributions due to the time-dependent correction, the thermody- 
namic factor and the shape of the profile. 



6. Effect on the ionization history and the CMB 
power spectra 

In this Section we now give the expected correction to the ioniza- 
tion history when including the processes discussed in this pa per. 
For this we modified the Recfast code dSeager et al.lll999l) . so 
that we can load the pre-computed change in the Sobolev escape 
probability studied here. 

In Fig. [17] we present the final curves for AP/P as obtained 
for the different processes discussed in this paper. In Fig. [18] 
we show the corresponding correction in the free electron frac- 



tion computed with the modified version of Recfast. One can 
clearly see that the dominant correction is due to the thermo- 
dynamic factor, resulting in AP/P ~ +10% at z ~ 1280 and 
AN e /N e ~ -1.9% at z ~ 1120. The next largest correction is 
due to the time-dependent aspects of the problem, leading to 
AP/P ~ -5.6% at z ~ 1120 and AP/P ~ +5.9% at z ~ 1490. 
The associated correction in the free electron fraction has a max- 
imum of AN e /N e ~ +1.2% at z ~ 1020. The smallest correc- 
tion is due to the exact shape of the effective line profile, re- 
sulting in AP/P ~ +2.0% at z ~ 1200 and AN e /N e ~ -0.4% at 
z ~ 1080. The total correction then is AP/P ~ +ll%atz ~ 1420 

and A/Ve/We 1.6% at z ~ 1190. This is an important speed 

up of hydrogen recombination, although at z ~ 1100 a large 
part of the correction due to f v alone is canceled by the time- 
dependent correction. At the maximum of the Thomson visibil- 
ity function z ~ 1100 we find AN e /N e ~ -1.3%, where about 
AN e /N e 0.4% is coming from the shape of the profile alone. 

For completeness we also show the correction that is due to 
the two-photon formulation alone, i.e. where we subtracted the 
time-dependent contribution from the total correction. It leads to 
AN e /N e ~-2.4%atz~ 1110. 

In Fig. [19] we finally show the changes in the CMB tem- 
perature and polarization power spectra coming from the total 
correction AN e /N e as given by the solid line in Fig. [18] In partic- 
ular the changes in the EE power spectrum are impressive, with 
peak to peak amplitude ~ 2% - 3% at I > 1500. One can also 
see that the modifications in the C/'s correspond more to a shift 
in the positions of the peaks rather than changes in the absolute 
amplitude. This is connected to the fact that the correction in 
the free electron fraction leads to a small shift in the maximum 
of the Thomson vi sibility, but does not affect the Silk damping 
length (ISilkl 119681) as much. It will be very important to take 
these changes into account in the analysis of future CMB data. 

7. Discussion and conclusions 

In this paper we gave a formulation of the Lyman a transfer 
equation which allows us to include the two-photon corrections 
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Fig. 19. Changes in the CMB temperature and polarization 
power spectra. We included effect of the time-dependent correc- 
tion, the thermodynamic correction factor and the profile correc- 
tion, simultaneously. We used the result for the 10 shell hydro- 
gen atom. 

for the 3s- Is and 3d- Is channels. We then solved this transfer 
equation and presented the results for the Lyman a distortion at 
different redshifts (Sect. [?}. From this we computed the effec- 
tive Lyman a escape probability and derived the corresponding 
corrections to the Sobolev escape probability (Sect. [5}. We used 
these results to compute the corrections in the free electron frac- 
tion during hydrogen recombination and the associated changes 
in the CMB power spectra (Sect.|SJ. Below we now shortly sum- 
marize the results of each of these Sections and also give a short 
discussion for future work and improvements. 

7. 1 . The resulting Lyman a spectral distortion 

In Sect. |4]we discussed the influence of the different processes 
on the Lyman or distortion at different redshifts. Including all 
the corrections considered here simultaneously one can conclude 
that at xd < the shape of the distortion is mainly determined 
by the time-dependence of the effective emission rate (cf. Fig. [5] 
and|6l right column). In the close vicinity of the resonance all 
sources of corrections under discussion here are important for 
the exact shape of the Lyman a spectral distortion at difference 
redshifts (cf. Fig. 2}. In particular inside the Doppler core the 
spectral distortion will be very close the value at the line center 
multiplied by l// v . On the blue wing the distortion is mainly 
determined by the shape of the line profile. 

With the method given here we can in principle also com- 
pute the Lyman a distortion, as it would be observable today. 
Since all the processes discussed here affect the exact shape of 
this distortion, one could in principle probe our understanding 
of the recombination dynamics by observing it. As mentioned 
earlier dChluba & Sunvaevll2008b|) . similarly one expects addi- 
tional changes in the detailed shape of the Lyman a distortion 
due to partial frequency redistribution and electron scattering. 
All these processes therefore not only affect the dynamics of cos- 
mological hydrogen recombination, but in principle should leave 
obs ervable traces in cosmolo gical recombination spectrum (e.g. 
see lSunvaev & Chlubal2007l and references therein) until today. 
Measuring the exact shape of the Lyman a distortion and the 
other recombinational lines would in principle allow us to reveal 



these differences, and therefore directly probe our understanding 
of the recombination problem. Also if something non-standard 
happens (e.g. due to early energy release by decaying or anni- 
hilating dark matter), this will affec t the exact shape of the cos - 
mological recombination radiation dChluba & Sunvaevll2008ah . 
Therefore, by observing the recombinational radiation one in 
principle can directly uncover potential unknowns in the cosmo- 
logical recombination problem, a task that may not be completed 
otherwise. 

7.2. Corrections to the escape probability 

We have discussed the corrections to the effective Lyman a es- 
cape probability, showing that the largest contribution is coming 
from the thermodynamic factor f v (see Fig.fTTI). The next largest 
correction is due to the time-dependent aspects of the recom- 
bination problem, where in the formulation g iven here the self- 
feedback problem dChluba & Sunvaevll2008bl) appearing at low 
redshifts when using the '1 + 1' photon picture disappears (cf. 
Fig. |T7T ). Furthermore the time-dependent correction partially 
cancels the correction for the thermodynamic factor at low red- 
shift (z < 1300 - 1400), leaving a significanly smaller change in 
the escape probability at z ~ 1 100. As we explained here, these 
two corrections can be also obtained in the normal '1 + 1' photon 
picture, but for the thermodynamic correction factor a natural in- 
terpretation can only be given within the two-photon picture (see 
Sect. I2.1.11 i. Here the crucial ingredient is that the spectrum in 
the vicinity of the second, low frequency photon y' is given by 
the CMB blackbody photon field, which then allows us to write 
the radiative transfer equation for the Lyman a photon as equa- 
tion for one photon only. A very similar formulation should be 
applicable in the case of expanding envelopes of planetary neb- 
ulae and stars, where the photon field in the vicinity of y' will be 
given by a weakly diluted blackbody spectrum. However, when 
the photon distribution in the vicinity of both photons involved 
in the absorption process are far from their equilibrium values 
the derived formulation does not work. 

We also showed that the correction coming from the exact 
shape of the line profile in the vicinity of the Lyman a reso- 
nance leads to the smallest separate correction under discussion 
here. Only this part of the correction is really related to quan- 
tum mechanical modification of the transfer problem. Still the 
final contribution related to this modification is significant at the 
required level of accuracy (cf. Fig.fTTb. 

7.3. Implications for the ionization history and the CMB 
power spectra and critical remarks 

The results for the changes in the free electron fraction and CMB 
power spectra are shown in Fig. [18] and [T9j respectively. The 
main effect on N s due to the processes discussed here is a net 
speed up of hydrogen recombination by AN e /N e ~ -1.3% at 

z ~ 1100. About AN e /N e 0.4% of this correction is coming 

from the shape of the line profile alone, while the rest is due to 
the thermodynamic factor and the time-dependent aspects of the 
recombination problem. Here we would like to emphasis again, 
that the latter two corrections can actually also be obtained in the 
standard '1 + 1' photon picture, when introducing the thermody- 
namic correction factor using the detailed balance principle. 

We would also like to mention that our results for the 
changes in the free elect ron fraction seem to be rather similar 
to those of lHiratal d2008l) for the contributions from high level 
two-photon decays alone. However, we obtained these with a 



24 



Chluba and Sunyaev: Ly a escape during cosmological hydrogen recombination 



completely independent method. For the future it will be very 
important to perform a more detailed comparison once all rele- 
vant additions are identified. 

Regarding the CMB power spectra, in particular the associ- 
ated changes in the EE power spectrum are impressive, reaching 
peak to peak amplitude ~ 2% - 3% at I > 1500 (see Fig. [19]>. It 
will be important to take these corrections into account for the 
analysis of future CMB data. 

However, it is also clear that several steps still have to be 
taken: (i) one still has to study more processes; and (ii) one has 
to device a sufficiently fast method to simultaneously incorpo- 
rate all the corrections discussed in the literature so far into the 
computations of the CMB power spectra, in order to be ready for 
real parameter estimations using the CMB data. 

Regarding the first point, for e xample , the co rrection due to 
Raman processes (as explained in Hirata (2008) mostly due to 
2s- Is Raman scattering) leads to an additional delay of recombi- 
nation at low redshifts with a maximum of AN e /N e — 1-0.94% at 
z ~ 920. We did not include this process here, but it certainly is 
very important for accurate predictions of the CMB power spec- 
tra and should be cross-validated in the future. It is clear that one 
should also include the effect of higher level two-photon decays 
(e.g. from the 4d-level), since they are expected to speed up hy- 
drogen recombination in addition, likely affecting the result ob- 
tained here by another- 10%-20% (i.e. ANJN e ~ 0.1-0.3%). 
And finally, the effects of partial frequency redistribution, line 
recoil and electron scattering should be studied. Here in partic- 
ular the effect of line recoil will be important, leading to a sys- 
tematic drift of photons towards lower frequencies which again 
accelera tes hydrogen recombination by AN e /N e ~ —1.2% a t 
z ~ 900 dGrachev & Dubrovichl2008l:IChluba & Sunyaevl2009l) 

Re garding the second point, fairly recently iFendt et ail 
d2008l) proposed a new approach called RicJn, which uses 
multi-dimensional polynomial regression to accurately represent 
the dependence of the free electron fraction on redshift and the 
cosmological parameters. Instead of running the full (slow) re- 
combination code, one trains Rico with a set of models, so that 
the interpolation between them will be very fast afterwards. This 
approach should allow to propagate all the corrections in the ion- 
ization history that are included in the full recombination code 
to the C MB power spectra , without using an y fudging like in 
Recfast dSeager et all 1 9991: IWong et al.ll2008l) . In the future, we 
plan to provide an updated training set for Rico, including the 
corrections discussed here. This should also make it easier for 
other groups to cross-validate our results and will allow us to 
focus the effort on the physics rather than on approximating it. 



Appendix A: Different channels for the emission 
and death of Lyman a photon involving a 
sequence of two transitions 

If we restrict ourselves to the main channels that can lead to the 
emission or absorption of photons in the vicinity of the Lyman 
a resonance and involve two photons then we are left with: (i) 
ns — > Is two-photon decay and absorption; (ii) nd — » Is two- 
photon decay and absorption; (iii) c — > Is two-photon recombi- 
nation and photoionization. 

The problem is now to compute the emission and absorp- 
tion profiles connected with these processes and describe their 
relative contributions or probabilities at a given frequency. This 
in principle can be done for all possible routes. However, here 



we will focus on formulating the problem for the 3s — > Is and 
3d — > Is two-photon channels, not altering all other channels. 
This is a reasonable first approximation, since as we show here 
it is already clear from the ' 1 + 1 ' photon picture that the main 
contribution to the death of photons is coming from the 3d chan- 
nel (see Fig. [TJ. We only also add the 3s channel, since it is 
conceptually very similar. 



A. 1. The 3s-1s two-photon channel 

In this section we would like to derive the rate equation that de- 
scribes the evolution of the population, N 3s , in the 3s level, but 
where we take the two-photon aspect of the 3s- Is transition into 
account. In the normal ' 1 + 1 ' photon picture this transition is 
described by the sequence 3s — » 2p — > Is and its inverse process 
Is — > 2p — > 3s. For the 3s rate equation it is therefore clear that 
the two-photon correction due to the 3s- Is channel should affect 
part of the 3s — » 2p and 2p — > 3s rate. Once this part is known 
one can in principle replace it in the rate equations using a more 
complete two-photon description. 

Here we just give the formulation of this problem, also be- 
cause it allows to understand the most important aspects of the 
two-photon picture. However, to compute the corrections to the 
escape probability we directly use the (pre-computed) solution 
for the populations given in the standard approach, and only 
solve for the presumably small correction in the evolution of the 
photon field around the Lyman a resonance. Solving the com- 
plete set of modified rate equations simultaneously will be the 
final goal, for which one will need the results of the derivation 
presented here. We leave this problem for a future work. 



A.1.1. Isolating the different T + 1' photon routes 

In order to isolate the contribution from the 3s- Is two-photon 
channel we start by writing down the 3 s rate equation in the ' 1 + 
1' photon picture, including all possible ways for electrons in 
and out of the 3s-level 



cUV 3s d/V 3s 



di 



dt 



3s2p 



+ *3s 



(A.l) 



Here we have directly separated the part due to the Balmer a 
transition 



d/V 3s 



dt 



3s2p 



- ^2p3s - ^3s2p N 3s 



83s 



= N 2 pA 3s2p — «pl(V32) - N 3s A 3s2p [1 + «pi(v 32 )] 
£2p 

(A.2) 

which below we want to discuss in more detail. Here iV; and g, 
denote the population and statistical weight of level i, « p i(v) is 
the CMB blackbody occupation number, and A,j and v, ; - are the 
transition rate and frequency between level ; and j. All the other 
possible channels in and out of the 3s-level lead to the terms 



^3s - ^3sc 



+ y. — A;3s «pl(Vi3 s ) 
^ 83: 



Rl = N e N p R c3s + J] NiA i3s [1 + « p i(v,3s)], 



(A.3) 
(A.4) 



i>3s 



http://cosmos.astro.uiuc.edu/rico/ 



where N e and N p are the free electron and proton number den- 
sities, and R C i and Rj C are the recombination and photoioniza- 
tion rates to the level ;. Note that, since at frequencies below 
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the Lyman a line the spectral distortions du ring recombination 
are small (e.g. cf. IChluba & Suny aev 2006aj), above we simply 
replaced «(v) — » « p i(v) everywhere. 

Now we are interested in refining the term connected with 
the Balmer a channel. Since any two-photon or multi-photon 
process only leads to tiny corrections in the total decay rates, it 
is possible to use the one photon rates to compute the relative 
contributions of different transition sequences. It is clear that the 
term N^A^~ = N 3s A 3s 2 V [1 + « p i(v32)] describes the total flow 
of electrons in the direction of the 2p-state. Once the electron 
reached there it can go back to the 3s level with the probability 



P2p3s - 



^3s2p fr «pl(>32) 

o-P 

A* +R- 

2p 1 s 2p 



(A.5) 



Here A 2pls = A 2p i s [l + « p i(v2i)] is the stimulated Lyman a emis- 
sion rate in the ambient CMB blackbody field, and R 2p is the total 
rate at which electrons can leave the 2p-level, but excluding the 
Lyman a line. It is given by 

R 2p = R 2 P c + y — Aa v « P i(v /2p ). (A.6) 

.?2 P 



/>2 P 



Similarly, the electron can take the route 2p 
ability 



P2pls 



^2 P ls 



^2 P ls 



R 



= l ? 

— Penr 



Is with the prob- 



(A.7) 



2p 



or any of the other channels (2p — * «s (« > 2), 2p — > nd or 
2p — > c) with probability p2pu - 1 ~ PipH- Note that here we 
have neglected the deviations from a blackbody in the stimulated 
Lyman a emission rate, which in any case is tiny. We just wanted 
to maintain the structure of the equations. 

Then it is clear that the term A^s A^ s2 can be interpreted as 



^3s A; s2p = N 3s A* 3s2p x [p 2pls + p 2p nl 



(A.8) 



because the total flow of electrons should split up into those elec- 
trons that go to the ls-level (probability /?2 P is) and those that 
don't (probability p2pu)- 

From the physical point of view it is now clear that the par- 
tial flow connected with p2puN 3s Aj s2p should be interpreted as 
3s- Is two-photon emission in the '1 + 1' -photon picture, which 
we will later replace with the more proper two-photon formulae. 
The rest (p2 P nN 3s Aj s2 ) describes the contributions of all the 
other possible channels, e.g. also including the Balmer a scatter- 
ing transition 3s — > 2p — > 3s. We will continue to describe all 
these in the '1 + 1' photon picture. 

In order to understand the term connected with the total flow 
of electrons from the 2p-level towards the 3s-state, given by 
A^2 P A3 S 2 P — « P i(v32), we have to think about an electron that is 
added to the 2p-state. It will take the route 2p — * 3s with prob- 
ability /? 2p 3s as given by Eq. ( |A.5t . If we consider all possible 
routes into the 2p-state, and write the corresponding total rate as 
R 2p ~ R 2p + -^ls A 2p i s «(v 2 i), then one has the identity 

A^2 P A 3s 2 P — «pl(V32) = P2p3s [^2 P + Ms — A 2p l s «(v 2 i)] 
$2p (?ls 



where 



Rl p = N e NpR c2p + 



^W/Ai2p[l +«pl(v,2 P )]. 

i>2p 



(A.9a) 



(A.9b) 



Now one can write 

dAT 3s 

3s2p 



At 
where 

dN 3 



dt 



i+i 



/?2p3s ^2 P ~ P2 P ls A^3s A5 s2p + 



82p 



dN 3s 



di 



3s2pls 



(A. 10a) 



= P2p3sMs A 2p ls n(V2l) - P2 P lsA^3s A 3s2p . 

3s2pls £ls 

(A. 10b) 



The first two terms in Eq. (lA.lOal ) describe the partial flow of 
electrons towards the 3s-state, but where it is certain that the 
electron did not pass through the Lyman a line before. The last 
term is the fractional contribution of the 3s <-> 2p «-> ls-channel 
in the 3s rate equation, but described in the '1+1' photon picture. 
This is the term which in the end we will want to replace with 
the two-photon formulae. 
If we now identify 

A3s2p A2pls 



^3sls ~ 



A* +R~ 

2pls 2p 



(A.ll) 



then we can finally rewrite Eq. dA.lObt as 



diV 3s 



dt 



= A lli s N u— «(v2i) n pi (v 32 ) 

3s2pls ' ' g\s 

-A^ ls AT 3s [l +« pl (v 2 i)][l +« P i(v 32 )]. (A. 12) 

Note that in vacuum one would have A,'. = A3 S 2 P , as it should 
be, since the electron will only have one way to leave the 2p- 
state. However, within an intense CMB background field, also 
the other channels will become active (e.g. 3s <-> 2p <-> c), so 
that part of the 3s «-> 2p flow will go through them. This reduces 
the effective decay rate A 



2y 



A.1.2. Replacing the 3s-1s channel in the '1 + 1' photon 
formulation with the two-photon expression 



We now want to replace the part due to 



with the more 



it I3s2pls 

proper two-photon terms. For this we have to ask the question 
how the 3s- Is two-photon term actually looks like. If one con- 
siders an electron that is initially in the 3s-state, then one can use 
the vacuum 3s- Is two-photon decay profile in order to derive the 
emission profile needed to describe the injection of Lyman a 
photons for the escape problem. Simple formulae for the nec- 
essary va cuum two-photon decay pr ofiles can be found in the 
literature ( IChluba & Sunvaevll2008ch . We shall normalize these 

profiles like dvdQ, = 1, where (f? y (v) already includes 

the motion of the atoms in the same way as for the normal 
Lorentzian lin es, usually leading to the Voigt-profiles (e.g. see 
lMihalasll978l) . For convenience we chose the Lyman a Doppler- 
width, Avd, in the normalization. 

With this the net change of the number density of electrons 
in the 3s level via the 3s- Is two-photon channel is given by 



d^V 3s 



dt 



2y 



3sls 



z(v3i - v) dv 



■4^3 



s J <p] y s (v)n(v) n( 

J ¥i£(v)[l + n(v)][l + n(v 3 i - v)] dv 



(A. 13a) 
A^ ls A^i s n(v2i)n p i(v32) 

«pi(v 2 i)][l +« P i(v 3 2)], (A. 13b) 
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where A?J. is the effective 3s-ls two-photon decay rate, which 
in vacuum i 



A3 S 2p, but within an ambient black- 
body radiation field should take the value given by Eq. ( IA. 1 Q . 
Furthermore, (p^ = <^(v)/Avd denotes the 3s- Is two-photon 
decay profile, and V31 is the 3s-ls transition frequency. 

For the approximation Eq. JA. 13bb three comments should 
be made: first we have assumed that the main contributions to 
the integrals over the two-photon line profiles are coming from 
the poles close to v ~ V21 and v ~ V32. Second, we have used the 
fact that the CMB spectral distortion around the Balmer a line 
are tiny. Also the stimulated term in the vicinty of the Lyman 
a resonance is completely negligable, so that we just can use 
1 +n(v2i) ~ 1 +n p i(v2i) instead, without changing anything. And 
finally, we assumed that only for the ls-3s two-photon absorp- 
tion rate the deviations of the CMB spectrum from a blackbody 
in the vicinity of the Lyman a resonance will matter. 

The result Eq. ( IA. 13bb is identical with the term given by the 
'1 + 1' photon picture, Eq. dA.12| ). This is not surprising, since 
with the above approximations we have simply turned from the 
two-photon to the '1 + 1' photon picture. In order to include the 
effect of two-photon transitions in to the rate equation of the 3s- 

level, we should therefore replace ^r^Lt, , with the full integral 

r at I3s2pls ° 

given by Eq. dA.13at . 



A.1 .3. Term in the Lyman a radiative transfer equation 



In order to use the integral ( |A.13at in the computations of the 
ionization history, we also have to give the solution of the CMB 
spectral distortion in the vicinity of the Lyman a resonance. We 
therefore have to explicitly write the 3s- Is two-photon emission 
and absorption terms for the evolution of the photon field and 
solve the corresponding transfer equation. In particular we want 
to bring the transfer equation into the form Eq. @. 
From Eq. (]A.13at it directly follows 



I ^1 
c dt 



2y 



3sls 



2A 



An 



A^(v)[l+«(v)][l + rt(v 3 i-v)] 



2A 2y 

-M*N ls <pll(v) «(v) «(v 3] - v). (A. 14) 



Here the factor of 2 is due to the fact that per electron two pho- 
tons are involved, and the factor of An converts to per steradian, 
having in mind that the medium is isotropic. Furthermore in this 
form it is assumed that every two-photon interaction in the 3s- Is 
channel leads to a complete redistribution of the photons over 
the whole two-photon profile. This also means that we have not 
distinguished two-photon emission and absorption from two- 
photon scattering events. However, this should be a very good 
approximation, since the scattering event involves two photons. 
This means that the total energy of the incoming photons will be 
split up such that in most cases the scattered photons will have 
energy v ~ V21 and V ~ V32. Note that this does not imply that we 
are using a complete redistribution approximation for the Lyman 
a resonance scattering itself, since only about ~ 10~ 4 - 10~ 3 of 
all interactio ns will lead to the 3s- and 3d-state via two-photon 
interactions (IChluba & Sunvaevll2008bT) . 



15 Here the approximate sign is due to the fact that the rate coefficient 
in the two-photon formulation should contain some small (~ 10~ 6 - 
10~ 5 ) quantum mechanical correction to the one photon rate. This will 
not lead to any significant correction in the escape probability. 



Neglecting the deviations from the blackbody spectrum in 
the emission term and comparing with Eq. <j4j we can identify 



»3s«l 



s(v) = 2 



v>v 31 /2 
i 



1 +n p i(v) 1 +«pi(v 3 i - v) 



l+«pl(v 2 i) l+«pl(V32) 
1 + «pi(v 3 i - V) 



1 + «pl(V3 2 ) 



(A. 15a) 



(A. 15b) 





A* 


ly 


/i 2pls 


Pern 


A 2pl S +*2p 


2p 





(A. 15c) 
(A.15d) 

This result shows that the effective profile of the process as ex- 
pected is given by the two-photon profile for the 3s- Is channel 
including the induced terms relative to the values at the Lyman 
and Balmer a resonance. Also the emission probability is ex- 
actly the Lyman a emission probability including the induced 
emission for the central frequency of the Lyman a line. And the 
last term simply represents the number density of 3s-electrons 
that reach the 2p-state per second in the '1 + 1' photon picture, 
where again the stimulated emission due to CMB photons close 
to the Balmer a frequency was included. 

Note that 03 SM i s (v) is no longer normalized to unity. In vac- 
uum one would find J^°° ^"^^ dv d£2 = 2, while within the 

CMB blackbody field £° dvdQ > 2. However the rel- 

ative correction to the overall normalization of the profile due 
to stimulated emission is of the order of ~ 1(T 9 - 1(T 8 (see 
Sect. 15.2. U . Also one should mention that due to the symme- 
try of the profile around v = V31/2, by restricting the range of 
integration to V31/2 < v < V31 one can avoid counting both the 
Lyman a and Balmer a photons. We will use this fact to simplify 
the numerical integration (see Sect. I2.3.T1 . 

If we now look at the absorption term in Eq. ( IA. 14t . using 
the definitions ( lA.15t we can directly write 



I ^1 
c ~dt 
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ls(v) ,2y . 

A 3s ' ls A^i s n(v)« p i(v 3 i - v) 



An Av D 

1 + «pl(V2l) 



1 + n p i(v 32 ) 



v) 



(A. 16) 



1 + « p i(v) 1 + n p i(v 3 i 

Here we have already assumed that the important part for our 
problem is the region v > V31/2. This implies that V31 - v < 
V31/2, so that the deviations from the CMB blackbody can be 
neglected leading to n(v)n p i(v3i - v) instead of «(v)n(v3i - v). 

Since according to Eq. CO]) and ( 1A.1U A] y sls = P2 P 3s^^j, 

and because [1 + n p i(v)]/n p i(v) = e hv ^ kT T, we finally find 



I ^1 

c dt 



2y 



_ 03 s ^l s (y)^2pls.P2p3s 
3sls,abs 47rAv D g3s/g2p 



l! 1 + «pl(v 2 l) 
1 + Hpl(v) 



"(v)Nls 

(A. 17a) 



r V3s^l s (v)g2p 

47rAv D g ls 



A2pisP2p3s e tT ? n{v)N u . (A.17b) 



Note that g3 S /g 2p = giJg2 P - 

With the Einstein relations it is then easy to show that 

j^- A2 P is n(v) = I1V21 Bu^rNy, so that we directly verify the 
thermodynamic correction factoi0 /v = -7- e h ^ V2 ^^ kTy and find 



16 From Eq. dA.17ab bycomparing with Eq. 01 we can see that with 
the choice of coefficients and variables we more rigorously infer f v = 
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^isjs = p 2p3s . We therefore have confirmed the completeness of 
the form of the Eq. © for the 3s- Is two-photon channel. 



A.2. The 3d-1s two-photon channel 

After going through the argument for the 3s- Is channel it is easy 
to do the same for the 3d- 1 s channel. For the rate equation analog 
to Eq. (|A. It . jA3\ , and dA.lOl ) one has 



dN M _ djV 3d 
dt dt 



3d2p 



+ Rt A - R M N M 



R M = R3dc + V, — A,3d "plO,3d) 

^3d = N e N p R cM + ^ N t A Bs [1 + n p i(v,3d)] 



where the Balmer a channel is defined by 
dN 3d 



dl 



3d2p 



1 + 1 

3d2pls 



P2 P 3d 

dN 3d 1+1 



(A. 18a) 
(A. 18b) 
(A.18c) 

(A.18d) 
(A.18e) 



dt 



3d2pls 



dN 3d 

P2p3d^2p - P2pls N 3d ^3d2p + _ ^7~ 
A 3d 2p fr p «pl( v 32) 

A* + R~ 

2pls 2p 
A 3dls ^ls— «( y 2l) «pl(K32) 



^dls^dfl +«pl(V2l)][l +«pl(^2)]. (A.18f) 



As before one should now replace Eq. ( IA. 1 8fb with 



dl 



= A 3dis— Ni* f ^(v)«(v)«(v3i-v)dv 



A3d2p/l2pl s 

A* + R~ 

2pls 2p 



(A. 19a) 
(A. 19b) 



if one is interested in the 3d- Is two-photon correction to the 3d- 
rate equation. Note that here the ratio of the statistical weights is 
not unity like in the case of the 3s- Is channel. 

The terms for the transfer equation can also be cast into the 
form <j4j where the important coefficients are given by 



03d«is(v) - 2(f> 3 ' A (v)- 



2^(v) 



l+«pl(v 2 l) l+« p l(V32) 
1 + « p l(V31 - V) 

1 + n p i(v 32 ) 



*S' + =^d2p^3d 
P d - /?2p3d- 



A.3. The other channels 



(A.20a) 
(A.20b) 

(A.20c) 
(A.20d) 



For the other channels in and out of the 2p-state we can also 
derive the corresponding partial rates in a similar way as for the 
3s and 3d state. However, since the main correction is expected 



± e -vf- = % S<21> « % Jtv-nMTy for hv » kT a condi- 

v 2 i+npiM i « P i(>') v 2 y 

tion that is fulfilled during cosmological recombination in the vicinity 
( v <i V31 /2) of the Lyman a resonance. 



to come from the 3s and 3d two-photon channels for these we 
will simply use the '1 + 1' photon picture. This means that we 
will not replace the corresponding ; <-> 2p rates with the two- 
photon description. For all the ns and nd-states with n > 3 the 
rate equations therefore will be similar to Eq. ( IA. 1 8t . and for 
electrons in the continuum one will have 



djV e 
dt 



dN e 
~~dt 



+ y [RicNt - N e N p R ci \ (A.21a) 

c2 p ;>2 



with 



dN, 
~d7 



C'2p 

P2pc 

dN e 1 + 1 



+ - dJV e 
= Pipe R 2 p _ /^pls N e N v R C 2 P + — j^- 

^2pc 



I + l 



c2pls 



A* + R~ 

2p 1 s 2p 



(A.21b) 
(A.21c) 



di 



82p 



c2pls 



Nis — A 2p i s «(v 2 i) pipe - Pipu N e N p R c2p . 

(A.21d) 



Still there is a small difference to the normal rate equations. In 
the formulation given above the population of the 2p-state has 
vanished from all the rate equations, and in particular from those 
for the 3s and 3d-state. For the 3s and 3d-state physically this 
is expected, since in the two-photon picture on the way to the 
ls-level the electron is not really passing through the 2p-state. 
In the full two-photon picture the electron reaches the Is level 
via all intermediate p-states, including those in the continuum. 
For the other levels the above formulation would have also been 
obtained by simply replacing the solution of the 2p-state with the 
quasi-stationary value in the '1 + 1' photon approach. In this way 
one again has a closed system of rate equations, which avoids the 
difficulty in attaching a population to the 2p-state. 

Appendix B: Derivation of the thermodynamic 
factor using the '1 + 1' photon picture. 

As mentioned in the introduction and also earlier 
dChluba & Sunyaevl l2008bl) . in the normal '1 + 1' pho- 
ton approximation the term describing the emission and 
absorption of Lyman a photons in full thermodynamic equi- 
librium is not exactly conserving a blackbody spectrum 
at all frequencies. This can be directly seen from Eq. (|T), 
since in full thermodynamic equilibrium one should have 
(Pem^pF = {p l d 7 "V2iB l2 N ls N v J\ Using the defini- 
tions of the previous Section, in equilibrium one expects 

KT - W ^ - «^ M l2 , g 4£ 

and p l Jn = A21 [1 + « P i(v2i)] P d y IR 2v > so ^ at one fi n d s 



lr hv 2 i B 12 N H N V ^' 



pIIRL 



2y\ 1 + « p l(V2l) 



(B.l) 



Here Ny — « P i(v) defines the blackbody spectrum. With the 
identity [1 + « p i(v)]/« p i(v) = e hv ^ kTy one therefore has 

'p/hvn B n N u N v ) eq _ 



r, 17 J? 



2p 



1 + « P i(v) *!a 
■ e k 



V*j 1 + « P i(v 2 i) 



l//v (B.2) 



The in the last step we assumed that hv » kT y and hv 2 \ » kT y , 
so that the factors [1 + «] could be neglected, an approxima- 
tion that is certainly possible during cosmological recombina- 
tion. However it in fact is only due to the used definition of 
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Pem> m which we evaluated the stimulated emission factors at 
the Lyman a line center only and in addition directly assumed 
a blackbody spectrum. More consistently one should have used 



Pern 



Pem [1 + n v], where we have 



ly.* 

Pern 



12pls 



A 2p i s [l +(n v ) em ]+R 2v 



(B.3) 



Here (n v ) em denotes the average of the photon occupation num- 
ber across the Lyman a emission profile. With this definition in 
Eq. ( lB.2t one would have obtained the factor l// y directly, with- 
out any additional comment. 

We would like to mention another way to obtain the ther- 
modynamic factor in the absorption process, which just uses the 
term in the Lyman a rate equation. In the standard textbooks 
(e.g. lMihalaslll978l) one finds: 



<W 2p 



N u — A 2 pis«(v2i)-iV2pA 2pls [l+n(v 2 i)]. (B.4) 



2pls 



Here two comments should be made: (i) more rigorously one 
should replace n(v 2 i) with the average occupation number over 
the line profile, i.e. «(v 2 i) — > (n v ) and (ii) in addition one should 
distinguish between the emission and absorption process, imply- 
ing that [1 + «(v 2 i)] — > [1 + (n v ) em ], while in the first term one 
has «(v 2 i) — > (n v ) a b s . In full equilibrium one should then find 
(n u ^A 2p h <ny> abs ) eq = (N 2p A 2pls [l +<«y> em ]) eq , so that 



f N ls g ±A 2?ls <« v ) abs ) 



,hv 2 i lkT y 



^ 2p A 2p i s [l + <n v )em] 

It is now easy to show that 



■(*). 



= 1 



(B.5) 



1+H 
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hv 2 ]/kT y 



abs 



My 2l -v\lkT y 



(B.6) 
(1+4) 

\ I e 



Since according to Eq. dB.5b (nf) e h ^-\i kT y 

\ / abs 

one should conclude that <n v ) abs = (e h ^ v ~ V21 ^ kTy n v ) , where 
n v now is an arbitrary photon occupation number. In terms of 

N v = ^n y one therefore has {N y ) llbs = <n v } a b s = {f v N v ) em . 
With this we obtained the thermodynamic correction factor f v in 
the absorption profile for N v - ~r n v , since (M<)abs = (fvN v ) em 
automatically implies a b s = f v <pem- 

Appendix C: Computation of two-photon profiles 

We com pute the two-photon deca y profiles according to the 
work of (IChluba & Sunvaevll2008bl) . There in particular the in- 
finite sum over intermediate states was split up into those states 
with principal quantum numbers n > «, and n < n,, where «,• is 
the initial state principal quantum number. This makes the sum 
over the resonances (in the case of 3s and 3d only one) finite and 
allows to give fitting formulae for the remaining contribution to 
the total matrix element coming from the infinite sum. This pro- 
cedure is very convenient for numerical evaluations. 

Here we would like to mention that the two-photon decay 
profiles behave like <p x v(y,i s ~ v) in the limits v — > or v — > 
v,-i s . This is because in this limit the main term in the infinte sum 
is coming from the Matrix element n,s/d — > n,p, which in the 



non-relativistic formulation has zero transition frequency. This 
implies that for v — > or v — > v,i s 

M G ntk y (1 - y) 1 3 2 n f* = ^ J («? - 1|4). (C. 1) 

Here 1 14 means 1 for the s- states or 4 for the d-states . Inserting 
numbers (for definitions see lChluba & Sunya ev 2008c) one finds 



*'' ls « 1-0598 xl0 4 y 



fa - l) 2n/ n? - 1 , 

v ~ s (C2a) 

*^ ls * 4.2393 x 10 3 y(l -y) ( "'' ~ 1^"' ^-^ s" 1 . (C.2b) 



(fU + l) 2n > Hi 



For the 3s and 3d profiles one therefore has 



^ ls (y)^441.6y(l-y)s- 1 
$Li.O0* 110.4); (l-^s" 1 



(C.3a) 
(C.3b) 



We will use these simple formulae to compute the two-photon 
spectra at < y < 0.001 and 0.999 < y < 1. 

The most important consequence of this limiting behavior 
with frequency is that due to stimulated emission in the ambi- 
ent CMB blackbody radiation field the two-photon profiles no 
longer vanish at y ~ and y ~ 1, since n p i(v) ~ 1/y for v ~ 0. 
For the 2s- Is two-photon proces s this behavior was also seen 
earlier (IChluba & Sunyaevl2006bl) . In the case of 3s and 3d two- 
photon decays this enhances the emission of photons close to 
the Lyman /3 resonance (cf. Fig. [2]). However, we find that the 
corrections due to stimulated two-photon emission are not im- 
portant for the cosmological recombination problem. 

Appendix D: Small corrections due to the motion of 
the atom 

To account for the motion of the atoms in the computations of 
the emission pr ofile one has t o compute the following integral 
(see Sect. 9.2 in lMihalaslll978l) 



0m(v) 



-f 



My'0))i 



dt 



(D.l) 



over the rest frame emission profile, r (v), which for a given 
frequency v due to the Doppler effect has to be evaluated at 
v'(f) = v[l - —t], where £ = c ^2kT/mu c 2 . The exponential 
factor arises from the Maxwell-Boltzmann velocity distribution 
for the neutral hydrogen atoms. 

For the Voigt-profile one normally uses the approximation 
v'(f) a; v - V21 — t, so that the emission profile can be written in 
terms of the Voigt-function 



H(x D , a) 



'dt 



(x D - t) 2 + a 2 



(D.2) 



for which simple app roximation in t erms of the Dawson-integral 
exist (see Sect. 9.2 in lMihalasll978h . Here a is the normal Voigt- 
parameter, and xu is the frequency distance from the line center 
in Doppler units of the Lyman a resonance. 

However, due to the steepness of the Lorentzian close to 
the line center this approximation leads to a small inaccuracy 
(A<f>/</> ~ 10~ 4 - 10~ 3 ), which actually is not necessary. To avoid 
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1 1 

Voigt <--> Voigt 

- 









Fig. D.l. Small correction to the normal Voigt-profile. We show 
the relative difference of the standard Voigt profile <p cc H(x D , a) 
in comparison with the modified expression based on H(xv, a) 
atz= 1100. 



this one should simply replace x D and a in the Voigt-integral 
with xd = xv/[l + xd^-] and a — a/[l + Xr>-^f], and in addi- 
tion multiply H with 1/[1 + xu—], yielding 



H(xd, a) = 



H(x D , a) 



II 



■x D ^]' 

U V 2 1 J 



As one can see that then the Voigt profile, <py 
effectively behaves as 



(D.3) 



H(xd, a) I y/x, 



/(xd) 



1 e~*v 
[1+^] 
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1 Av °n . o 2x 
1 -x D (1 + 2x D ) 



(D.4) 



close to the line center. At xd ~ ±1 this implies a relative cor- 



rection of A<p/(p 



=6 x 10~ 3 at z ~ 1100. However, 



at xd ~ 3 this corrections is expected to reach the 0.1% level. In 
Fig. ID. li the frequency dependence of this correction is shown in 
more detail, confirming these statements. Note that as expected 
the behavior in distant wings Qxo\ » 1) is not changed. 



D. 1 . Two-photon profiles for moving atoms 

As mentioned above, to include the motion of the atoms in 
the computation of the two-photon profiles one can in princi- 
ple directly use the expression (ID. lb . However, the computa- 
tion of this integral is rather time-consuming, and in the very 
distant wings also is not necessarjQ- Therefore we use a dif- 
ferent approach in which we utilize the fact that the relative 
difference, A<^(v)/0a(v), of the restframe two photon profile, 
<p ri (v), and the normal Lorentzian, <f>\(v), varies very slowly 



17 There the motion of the atoms can be neglected since the two- 
photon profiles normally vary very slowly with frequency. Although 
this is not entirely true for the ns two-photon profiles close to the fre- 
quencies where 2y vanishes, one expects a negligible additional cor- 
rection due to this approximation. 



3s <— > Lorentz 

3d <— > Lorentz 

3s <-> Voigt 

3d <-> Voigt 

analytic approximations 




Fig. D.2. Relative difference of the 3s and 3d two-photon pro- 
files with respect to the Lorentzian (thin line) and the Voigt pro- 
file (thick lines). In the first case the motion of the atoms was 
neglected, while for the comparison with the Voigt profile it was 
included. The boxes show the analytic approximation for the cor- 
rection including the motion. 

Table D.l. Parameters for the 3s and 3d two photon profiles 
around the line center. These were determined in the range 
-10 < xd < 10. For explanation see Appendix ID1 







Pi 


3s 


3.73335 x KT" 


+6.5325 x 1(T 4 


3d 


7.45559 x 1(T 7 


-1.5129 x 1(T 4 



with frequency, as compared to 0a (v). Therefore, to lowest or- 
der for moving atoms the two-photon profile will be given by 

<fyxD,a) * <Mx D ,a) [ l + A^J(v)/0a(v)]. 

However, close to the line center this approximation due 
to the steepness of the Lorentzian again becomes inaccurate at 
the percent level. The lowest order correction can be computed 
approximating the rest frame two-photon profiles close to the 
Lyman a line center by 



#00 



y(v)[l+ai +PiX D ]. 



(D.5) 



Here is the The coefficients a, and for the 3s and 3d level are 
given in Table ID. f| 

Us ing the same method as described in Sect. 9.2 of lMihalasI 
dl978l) it is then easy to show that for |x D | < 1000 



/,-}' 



(x D ,a)x(f> v (x D ,a) 



A0j(v) pi d rxD H(x D ,a) 



1 + 



<f> A (v) 2 H(xD,a) 



(D.6) 



Close to the line-center we will use this expression for the two- 
photon profiles. However, far away from the Lyman a resonance 
we will neglect the effect due to the motion of the atoms, and 
simply use the rest frame two-photon profiles. 

However, we find that the correction in connection with the 
motion of the atom are not important for the cosmological re- 
combination problem. 
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